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ABSTRACT 


The  methods  for  generating  and  testing  uniform 
pseudo-random  numbers  are  reviewed.  Some  techniques 
for  transforming  uniform  random  deviates  to  non-uniform 
random  deviates  are  described,  principally  for  the 
normal  distribution. 

Applications  show  some  uses  to  which  sequences  of 
pseudo-random  numbers  can  be  adapted.  A  simulation  is 
introduced  which  examines  probability  estimates  for  a 
binomial  distribution. 

The  Appendix  shows  results  for  some  tests  which 
indicate  that  multiplicat ive-congruent ial  pseudo-random 
sequences  are  a  satisfactory  source  of  random  numbers. 
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CHAPTER  I 


INTRODUCTION 


Many  investigations  involving  applied  mathematics, 
statistics,  or  physics  make  use  of  Monte  Carlo  methods. 

For  example,  in  applied  mathematics,  we  may  wish  to  obtain 
an  approximation  to  the  area  under  a  curve,  i.e.,  calculate 
J1  f(x)dx  .  We  could  select  a  pair  of  co-ordinates  (x,y) 


0 

at  random  and  test  the  value  of  f(x)  against  the  value  of 
y  .  If  f(x)  <_  y  ,  then  we  would  accept  the  point  as  falling 
on  or  below  the  curve;  if  f(x)  >  y  ,  we  would  reject  the 
point.  Finally,  we  would  compare  the  ratio  of  the  number 
of  accepted  points  to  the  total  number  of  points  selected 
at  random.  The  ratio  would  represent  an  approximation  to 
the  area.  Such  a  method  would  only  give  a  very  rough 
approximation  to  the  area  and  would  be  used  mainly  in 
applications  where  a  solution  could  not  be  obtained 
analytically.  In  physics,  we  could  use  Monte  Carlo 
techniques  to  examine  the  behaviour  of  the  diffusion  of  a 
gas,  simulate  random  collisions  of  a  molecule,  or  examine 
the  shielding  effects  of  a  substance  such  as  water.  One 
very  useful  example  is  the  simulation  of  the  path  of  a 


particle  with  Brownian  movement.  In  statistics,  we  might 
wish  to  simulate  a  game  of  chance,  such  as  coin  tossing 
or  a  queueing  problem.  In  all  the  preceding  examples, 


, 
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Monte  Carlo  techniques  would  be  used.  Such  techniques 
depend  upon  having  available  sequences  of  numbers  which 
appear  to  have  been  drawn  at  random  from  a  particular 
probability  distribution. 

Monte  Carlo  techniques  are  able  to  give  at  least 
approximate  answers  where  other  techniques  fail.  For 
example,  an  experiment  to  examine  absorption  of  X-rays 
would  be  very  difficult  to  control  and  to  obtain  any 
results  from.  However,  with  Monte  Carlo  methods,  we  would 
simply  generate  a  sequence  of  random  numbers  to  follow 
the  history  of  an  individual  ray,  perform  this  experiment 
for  a  sufficiently  large  number  of  trials  and  tabulate 
the  results. 

It  is  the  aim  of  this  thesis  to  examine  sequences 
of  random  numbers  generated  by  digital  computer  programs 
and  to  report  results  pertaining  to  some  special  sequences. 

Chapter  II  examines  the  methods  which  have  been 
used  in  the  past  and  current  methods  for  generating 
uniformly-distributed  sequences  of  random  numbers.  It 
shows  the  advantages  and  disadvantages  of  several  methods 
for  generating  random  sequences  on  computers  from  the 
point  of  view  of  statistics  and  the  speed  with  which  the 
sequences  are  generated. 

Chapter  III  surveys  different  statistical  tests  for 


# 
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examining  sequences  of  numbers  for  randomness.  Both 
methods  used  in  the  past  and  current  methods  are  examined. 
The  emphasis  is  placed  on  four  tests  devised  by  Kendall 
and  Babington-Smith . 

Chapter  IV  considers  the  problem  of  transforming 
a  uniformly  distributed  random  variable  into,  particularly, 
a  normally  distributed  random  variable.  Five  different 
transformations  are  illustrated  and  their  advantages  and 
disadvantages  are  discussed.  The  transformation  to  an 
exponentially  distributed  random  variable  and  a  Poisson 
variable  are  treated  briefly. 

Chapter  V  discusses  four  simple  applications,  three 
from  probability  theory,  and  one  from  mathematics. 
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CHAPTER  II 


A  SURVEY  OF  METHODS  FOR  GENERATING 
RANDOM  NUMBERS 


2 . 1  Introduction 

Student  (1908)  appears  to  have  been  the  first  to 
use  random  sampling  techniques  to  estimate  distribution 
functions.  His  method  of  choosing  a  random  number 
consisted  of  using  a  correlation  table  of  the  heights  and 
the  left  middle  finger  measurements  of  3000  criminals  from 
a  paper  by  MacDonell  (1901) „  The  digits  from  the  table 
were  written  on  3000  pieces  of  cardboard,  shuffled,  and 
were  drawn  at  random,  with  replacement,  four  at  a  time. 
This  method,  however,  proved  to  be  very  slow,  and  it  was 
very  difficult  to  determine  when  the  pieces  of  cardboard 
had  been  shuffled  well.  Karl  Pearson  suggested  to  Tippett 
(1925)  that  he  should  replace  the  entire  system  of  ticket 
sampling  by  a  table  of  random  numbers  ranging  from  0000 
to  9999.  Tippett's  table  (1927)  was  formed  by  taking 
90,000  digits  "at  random"  from  census  tables.  They  were 
grouped  in  fours  to  give  the  required  10,000  numbers. 

Karl  Pearson  has  shown  in  the  foreword  how  these  numbers 
can  be  converted  to  give  random  samples  from  a  non-uniform 
distribution.  Fisher  and  Yates  (1938)  produced  a  table 
of  15,000  random  sampling  numbers.  Their  table  was 
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compiled  from  among  the  fi f teenth-t o-nineteenth  digits  in 
certain  sections  of  Thompson’s  " Logarithmica  Brittanica". 
Fisher  and  Yates  applied  the  frequency  test  to  their 
table  of  15,000  digits  <>  They  found  an  excess  of  the 
digits  3,  6,  and  9  as  shown  in  Table  2.1.1.  If  the  digits 
were  uniformly  distributed,  then  the  expected  frequency 
of  each  digit  would  be  1500, 


Digit 

Frequency 

0 

1493 

1 

1441 

2 

1461 

3 

*  1552 

4 

1494 

5 

1454 

6 

*  1613 

7 

1491 

8 

1482 

9 

*  1519 

Table  2.1.1  - 

Observed  Frequencies 

Fisher  and  Yates  Digits 

In  order  to  correct  the  table,  they  removed  fifty  of  the 
sixes  "strictly  at  random"  and  replaced  them  by  one  of 
the  other  nine  digits  "strictly  at  random" 0  The  tables 
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then  satisfied  the  tests  applied  by  Fisher  and  Yates. 

2 . 2  Mechanical  Methods  for  Generating  Random  Numbers 

Kendall  and  Babington-Smith  (1938,1939a)  had 
considered  forming  random  sampling  numbers  from  mathematical 
tables  but  decided  against  it  in  favour  of  a  mechanical 
method,  There  were  known  to  be  non-random  properties  in 
sequences  of  digits  from  mathematical  tables,  such  as 
tables  of  logarithms.  As  an  example,  Kendall  and  Babington- 
Smith  cite  the  following  theorem  proven  by  Franel  (1917) » 

Theorem  2.2,1:  Consider  the  logarithms  to  base  10  of  the 
natural  numbers  from  1  onwards.  The  proportional  frequency 
of  any  digit  in  this  series  does  not  tend  to  a  limit. 

Such  a  series  of  numbers  will  contain  an  increasing  number 
of  runs  of  certain  digits. 

Kendall  and  Babington-Smith  were  the  first  to 
successfully  generate  random  numbers  by  means  of  a 
mechanical  device.  Previously,  mechanical  methods  had 
been  considered  untrustworthy.  However,  Kendall  and 
Babington-Smith ’ s  machine  was  designed  to  eliminate  the 
sources  of  bias  which  had  appeared  in  other  generators. 

Their  randomizing  machine  was  composed  of  a  disc  which 
was  rotated  by  an  electric  motor  at  a  very  rapid  rate  in 
a  darkened  room.  The  disc  was  divided  equally  into  ten 


* 

3- 


7 


sectors  on  which  appeared  the  digits  0  through  9, 
inclusive.  An  electric  spark  or  a  neon  lamp  illuminated 
the  disc  instantaneously  so  as  to  make  the  disc  appear 
stationary  for  a  moment.  When  the  flash  occurred,  a 
number  was  selected  from  the  disc  by  means  of  a  fixed 
pointer.  The  intervals  of  the  flashes  were  varied  by 
means  of  a  neon  lamp  in  parallel  with  a  condensor  in  an 
independent  electrical  circuit.  To  add  to  the  randomness, 
a  key  tapped  by  an  observer  broke  the  circuit 
intermitt ant ly  at  irregular  intervals.  The  table  of 
random  numbers  generated  by  this  mechanism  (see  Kendall 
and  Babingt on-Smith  (1939b))  satisfied  the  requirements 
for  random  numbers  for  approximately  ten  years. 

In  addition  to  constructing  a  table  of  random 
numbers,  Kendall  and  Babingt on-Smith  devised  a  series  of 
tests  for  sequences  of  random  numbers.  These  tests  will 
be  described  in  the  next  chapter. 

The  most  extensive  table  of  random  digits,  to  date, 
has  been  published  by  the  RAND  Corporation  (1955)-  The 
random  digits  were  produced  by  a  randomization  of  a  basic 
table  generated  by  an  electronic  roulette  wheel.  The 
process  required  a  frequency  pulse  source  providing,  on 
the  average,  about  100,000  pulses  per  second.  These 
pulses  were  gated  approximately  once  per  second  by  a 
constant  frequency  pulse.  Pulse  standardization  circuits 


■ 


passed  the  pulses  through  a  five-place  binary  counter.  In 
principle,  the  machine  was  a  thirty-two  place  roulette 
wheel  which  made  approximately  3000  revolutions  per  trial 
and  produced  one  number  per  secondo  A  binary-t o-declmal 
converter  was  used  to  convert  twenty  of  the  thirty-two 
digits  produced  and  only  the  final  digit  of  the  two-digit 
decimal  number  was  retained.  The  final  digit  was  then 
punched  on  a  card.  Production  of  the  digits  began  in 
April  of  194fo  After  500,000  digits  were  produced,  tests 
were  performed  on  them  with  satisfactory  results.  Later, 
after  continuous  running  of  the  generator  for  more  than  a 
month,  tests  showed  that  there  seemed  to  be  a  slight 
tendency  favouring  even  digits  more  than  odd  digits.  From 
this,  it  was  concluded  that  the  machine  had  probably  been 
running  down  during  the  month.  In  order  to  correct  the 
fault,  it  was  decided  to  randomize  the  digits  produced  for 
the  table.  Each  of  the  fifty  digits  punched  on  a  card  was 
added,  modulo  10,  to  the  corresponding  digit  of  the 
previous  card  and  the  result  punched  on  another  card. 

Tests  performed  on  the  derived  series  yielded  acceptable 
results  and  were  adopted  as  the  table  of  random  digits 
which  was  published. 

2.3  Random  Number  Generation  with  an  Automatic  Computer 


With  the  introduction  of  computers  during  the  late 
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1940’s,  the  use  of  tables  of  random  digits  became  almost 
impossible  and  was  certainly  undesirable.  The  tables 
required  large  amounts  of  space  for  storage,  and  access 
to  auxiliary  storage  devices  by  the  computer  was  slow. 

In  order  to  be  able  to  use  random  digits,  it  was  necessary 
to  derive  some  method  of  generating  them  within  the 
computer  itself.  Two  types  of  generators,  physical  and 
arithmetic,  have  been  suggested  for  internal  computer  use, 
and  will  be  discussed  in  the  remaining  sections  of  this 
chapter.  The  emphasis  will  be  on  arithmetic  generators. 

2  .  4  Physical  Random  Number  Generators 

Basically,  a  physical  random  number  generator 
consists  of  some  external  device  which  delivers  a  series 
of  random  pulses  to  the  computer.  The  computer,  in  turn, 
transforms  the  series  of  pulses  into  a  random  number. 

There  are  two  main  methods  of  obtaining  random  numbers 
by  means  of  a  physical  process.  The  first  method  is 
based  on  the  noise  of  electronic  tubes,  and  the  second 
on  the  radiation  of  radioactive  substances. 

The  basis  of  the  first  method  is  a  noise  generator. 
In  electronic  circuits,  there  is  inherent  fluctuating 
noise,  which,  with  appropriate  amplificaton ,  will  ensure 
a  fluctuation  in  the  output  voltage.  One  particular 
noise  generator  circuit  described  by  Shreider  (1964)  uses 
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a  gas-discharge  tube  and  a  magnet.  Noise  pulses  can  be 
obtained  directly  from  the  gas-discharge  tube  with  a 
suitable  orientation  of  the  magnet.  The  output  signals 
from  the  tube  are  applied  to  the  input  of  an  amplifier. 

Radioactive  sources  for  random  number  generation 
usually  consist  of  a  radiation  source  of  radioactive 
particles  plus  a  counter.  The  counter  registers  the 
number  of  radioactive  particles  released  in  an  interval 
of  time  At  .  If  the  number  of  particles  is  even,  the 
digit  zero  is  recorded;  if  the  number  of  particles  is 
odd,  the  digit  one  is  recorded  as  the  value  of  the  random 
digit . 

Pawlak  (1956)  has  used  a  random  number  generator 
based  on  the  random  state  of  an  electronic  circuit.  He 
lets  A  and  B  represent  two  possible  states  of  a 
flip-flop.  If  the  contact  is  switched  on,  the  flip-flop 
will  be  randomly  at  one  of  its  states  A  and  B  .  From 
this  base,  we  may  get  a  sequence  of  2k  random  elements. 
This  will  give  us  a  random  sequence  of  states  A  and  B 
which  are  statistically  independent .  One  of  the  sets 
of  positions  of  the  flip-flop  after  switching  may  be: 


ABBAAABBABBBAABAB . 
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One  quite  well  known  physical  random  number 
generator  is  ERNIE  (Electronic  Random  Number  Indicator 
Equipment)  described  by  Thompson  (1959).  The  machine  is 
used  in  the  British  General  Post  Office’s  Premium  Savings 
Bond  scheme.  In  the  scheme,  certain  bond  numbers  are 
eligible  to  take  part  in  a  draw,  and,  if  the  number 
generated  by  ERNIE  coincides  with  one  of  the  bond  numbers,  a 
prize  is  awarded  to  the  owner  of  that  bond  number.  In  the 
first  stage,  ERNIE  generates  and  prints  a  sample,  with 
replacement,  from  the  population  of  eligible  bond  numbers. 

In  the  second  stage,  reference  is  made  to  the  records, 
and,  if  a  bond  number  which  appears  has  already  been 
awarded  a  better  prize  than  the  current  one,  the  number 
is  deleted  from  the  list,  i.e.,  only  one  prize  can  be  won 
at  a  time.  Each  bond  number  contains  nine  digits  and 
therefore  ERNIE  must  generate  nine  random  digits.  In 
order  to  do  this,  the  electrical  noise  in  ten  neon  tubes 
is  used.  The  bond  numbers  produced  in  any  run  depend  on 
the  particular  noise  waveforms  produced  by  the  ten  neon 
tubes.  Figure  2.4.1  shows  the  parts  of  a  digit  generator. 
Each  time  the  noise  waveform  passes  upwardly  through  a 
certain  fixed  level,  the  pulse  generator  emits  a  short 
pulse  of  approximately  two  microseconds.  After  each  pulse 
is  generated,  there  is  a  delay  time  of  thirty  microseconds 
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before  another  pulse  can  be  generated  by  the  noise.  A 
train  of  pulses  is  thus  passed  through  an  n-position 
cyclic  counter,  each  pulse  advancing  the  counter  by  one. 
The  counter  is  examined  at  regular  intervals,  and  the 
digit  indicated  is  passed  to  the  output  as  one  of  a 
sequence  of  random  digits.  In  order  that  no  fault  should 
develop  during  the  selection  of  bond  numbers,  a 
redundancy  technique  has  been  incorporated  in  the 
generation  of  random  digits.  Ten  cyclic  counters  are 
provided  to  give  nine  random  digits.  The  counters  are 


Noise  Source 


Pulse  Generator 


Cyclic  Counter 


Inspection  Signals 


Sequence  of  Random  Digits 


Figure  2.4.1  -  ERNIE  Digit  Generator 
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connected  In  pairs  and  the  final  outputs  are  taken  as  the 
differences  of  each  pair  of  generators.  It  might  be 
pointed  out  that  statistical  tests  performed  on  ERNIE 
satisfied  the  criteria  for  randomness.  However,  the  use 
of  ERNIE  as  a  general  purpose  generator  is  very  restrictive 
since  the  speed  of  generation  is  extremely  slow.  For  the 
application  it  was  designed,  ERNIE  is  excellent  because 
relatively  few  bond  numbers  are  required  at  a  time. 


2 . 5  Arithmetic  Random  Number  Generators 

Arithmetic  generators  all  have  the  same  basic 
property  -  two  numbers  are  multiplied  together  and  some 
part  of  their  product  is  retained  as  a  random  sequence 
of  digits.  The  distinct  advantage  with  arithmetic 
generators  as  compared  to  the  physical  generators  is  that 
any  sequence  of  digits  produced  can  be  reproduced 
exactly,  if  required.  This  enables  one  to  correct  an 
error  in  a  program  and  use  the  same  data  to  check  the 
correction,  i.e.,  only  the  correction  and  not  a  different 
set  of  random  numbers  will  affect  the  result.  Since  all 
sequences  of  random  digits  produced  by  an  arithmetic 
method  can  be  predicted,  they  will  be  referred  to  as 
pseudo-random  numbers.  They  will  be  "random"  so  far  as 
satisfying  the  criteria  for  randomness  discussed  in  the 
next  chapter  is  concerned. 
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The  first  arithmetic  generator  was  proposed  in 
1946  by  von  Neumann  and  Metropolis  (see  Richtmeyer  ( 19 6 1 ) ) . 

This  was  the  middle-of-the-squares  process.  An  n-digit 

2 

number  Xq  is  squared  to  give  a  2n-digit  product  Xq  . 

The  middle  n  digits  are  retained  as  x1  and  the  process 
is  repeated.  As  an  example,  consider  Xq  =  1234  ,  then 
x^  =  01522756  and  x1  =  5227  .  x^  =  27321529  and 
therefore  x ^  =  3215  •  Tests  performed  on  the  numbers 
generated  by  the  middle-of-the~squares  method  produced 
unsatisfactory  results.  The  length  of  a  cycle  is  dependent 
on  the  starting  value  and  some  initial  values  can  lead  to 
short  cycles.  For  example,  with  xQ  =  3600  ,  the  sequence 
9600,  1600,  5600,  3600,  9600,...  is  produced.  Cycles  of 
this  type  are  difficult  to  detect  as  they  may  occur  at  any 
point  during  the  generation  of  the  sequences.  Another 
drawback  of  the  middle-of-the-squares  is  that  certain 
values  can  lead  to  a  degenerate  cycle,  i.e.,  a  result 
where  the  length  of  the  cycle  is  one.  Apart  from  the 
preceding  problems,  the  middle-of-the-squares  generator 
does  not  give  a  uniform  distribution  of  digits  regardless 
of  the  starting  value. 

An  improvement  to  the  middle-of-the-squares 
generator  is  the  mid-product  method  (see  Forsythe  (1951)). 
Initially,  two  n-digit  values  are  needed.  The  two  values 


f 
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are  multiplied  together  to  form  the  succeeding  number 
and  the  process  is  repeated.  For  example  xQ  =  1234  , 
x1  =  5678  give  the  series  0066,  3747  ,  2473,  2663,... 

The  advantage  of  the  mid-product  generator  as  compared 
to  the  middle-of-the-squares  generator  is  that  a  much 
longer  cycle  may  be  obtained.  Both  initial  values  must 
occur  consecutively  before  the  cycle  can  repeat. 

However,  this  method,  too,  has  values  which  can  lead  to 
a  degenerate  cycle  of  unit  length.  For  example,  with 
Xq  =  0100  ,  x^  =  0500  ,  the  series  is  0500  ,  2500  ,  2500  , 
2500 ,..  . 

Due  to  the  finite  word-length  of  a  computer,  any 
sequence  of  random  digits  must  eventually  repeat  itself. 
The  problem,  then,  is  to  make  the  cycle  as  long  as 
possible.  A  solution  was  proposed  by  Lehmer  (1949).  His 
method,  to  date,  has  been  the  most  successful  and  widely- 
used  for  generating  pseudo-random  numbers  on  a  computer. 
It  eliminates  any  possibility  of  a  sequence  degenerating 
into  a  small  loop. 

In  order  to  be  able  to  understand  Lehmer’ s  method, 
we  must  first  consider  the  definition  of  primitive  roots 
as  given,  for  example,  in  Abramowitz  and  Stegun  (1964). 

Definition  2.5.1:  The  integers  not  exceeding  and 


relatively  prime  to  a  fixed  integer  n  form  a  group;  the 


' 7  ° ;  > 


. 
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group  is  cyclic  if  and  only  if  n  =  2  or  4  or  n  is  of 

the  form  pk  or  2pk  ,  where  p  is  an  odd  prime.  Then 

g  is  a  primitive  root  of  n  if  it  generates  that  group, 

2  2  $  ( n ) 

i.e.,  if  g,  g  ,  g  ,...,  g  are  distinct  modulo  n  . 

There  are  $($(n))  primitive  roots  of  n  ,  where  $(n) 
is  the  number  of  integers  not  exceeding  and  relatively 
prime  to  n  (Euler's  <J>  function). 


The  results  for  selected  integers  are  shown  in  Table  2.5.1 


n  $(n)  $($(n))  Relatively  Prime 

Integers 

11  10  4  123456789  10 

9  6  2  124578 

18  6  2  1  5  7  11  13  17 

10  4  2  1379 

15  8  4  1  2  4  7  8  11  13  14 


Primitive  Roots 
of  n 

2678 
2  5 
5  11 
3  7 
NONE 


Table  2.5.1  -  Primitive  Roots  of  Selected  Integers 


For  the  derivation  of  a  method  for  finding  primitive 
roots  of  primes,  consult  Ore  (1948).  For  powers  of  2,  the 
primitive  roots  are  r  =  1  ,  and  r  =  3  (mod  4)  .  For 
modulo  8  and  higher  powers  of  2,  no  primitive  roots  exist 
because  all  odd  numbers  are  relatively  prime  to  the  moduli 
and  have  the  form 


a 


2n  +  1  . 


mnol  srlJ  jvsd  brB 
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Therefore 

a2  =  4n2  +  4n  +  1  =  4(n+l)n  +  1  . 

One  of  n  or  n  +  1  must  contain  a  factor  of  2,  so  that 

a2  =  1  (mod  8)  . 

But  $(8)  =  4  and  therefore  has  no  primitive  roots. 

Since 

a2  =  1  +  8t  , 

then 

a  ^  =  1  (mod  16)  , 
a^  =  1  (mod  32)  , 


and 


a2a"2  e  l  (mod  2a )  .  (2.5.1) 

Since  $(2a)  =  2a  ^  ,  (2.5.1)  shows  that  there  can  be  no 
primitive  roots  for  powers  of  2  higher  than  a  =  2  . 
Equation  (2.5.1)  also  implies  that  the  highest  exponent 

cx  a  “ *  2 

to  which  a  number  can  possibly  belong  (mod  2  )  is  2 


if  a  >  3  . 
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By  a  theorem  in  number  theory,  among  the  powers 
of  2,  only  2  and  4  have  primitive  roots.  For  all  higher 
powers  of  2a  ,  a  >_  3  s  every  odd  number  satisfies  the 
congruence 


a 


2^ 


1  (mod  2  ) 


(2.5.2) 


Lehmer's  method,  or  the  multiplicative-congruential 
method,  forms  a  sequence  of  pseudo-random  numbers 
according  to  the  formula 

xi+1  =  a  xi  (mod  M)  (2.5.3) 


for  given  a  and  subject  to  certain  restrictions 

which  must  be  observed  to  ensure  a  maximal  period.  These 
restrictions  may  be  summarized  in  the  following  theorem: 

Theorem  2.5.1:  The  sequence  defined  by  the  congruence 
relation  (2.5.3)  has  maximal  period  provided  that 
i)  Xq  is  relatively  prime  to  M  , 
ii)  a  is  a  primitive  root  for  pa  ,  if  pa  is 
a  factor  of  M  ,  with  p  odd  and  a  as 
large  as  possible. 

or,  with  p  =  2  and  a  =  1  or  2  . 


>- 
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lii)  a  belongs  to  2a  ^  if  2a  is  a  factor  of 
M  ,  with  a  >_  2  .  Moreover,  for  any  value 
of  M  ,  there  exist  values  of  a  satisfying 
these  conditions,  and,  finally,  the  maximal 
period  is  the  lowest  common  multiple  of  the 
periods,  (p-l)pa  ^  or  2a-^  with  respect 
to  the  prime  power  factors. 

It  is  relatively  easy  to  satisfy  condition  (i). 
Condition  (iii)  may  be  satisfied  if  a  =  ±3  (mod  8)  . 

To  show  that  this  is  true,  let  us  follow  an  example  with 
M  =  2P  .  In  this  case,  a  must  be  odd  to  be  relatively 
prime  to  2P  ,  i  .  e  .  ,  a  =  2m+l  .  Repeated  application 
of  formula  2.5°  3  gives 

xi+n  =  an  xi  (mod  2P)  (2.5.9) 

The  period  of  the  sequence  is  defined  to  be  the  smallest 

n  id 

integer  nQ  for  which  a  =  1  (mod  2F )  and  we  have 

p 

already  seen  that  n^  for  2P  is  2P”  .  The  general 

form  of  an  is 

an  =  (l+2m)n  =  1+n  2m+ ( ^ ) ( 2m) 2  +  (^ ) ( 2m) ^ ) ( 2m) V . . .  (2. 5. 5) 
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_  2 

For  n  =  2P  ,  (2.5*5)  becomes 


>P-2 


a  =  1+2F  2m+ 


p-20.„,2p  2(2p  2-1),„_.n2.2p  2(2p  2-1)(2p  2-2)/0_n3 


2! 


-( 2m)  +- 


3 


■(  2m) 


+  higher  powers  of  at  least  21 


)P-2 


e  1+2P  2m-2p  ^m2  (mod  2P) 


e  1+2P  ^m(l-m)  (mod  2P)  . 


(2.5*6) 


If  m  is  odd,  let  m  =  2b+l  .  Then  (2.5*6)  becomes 


a 


p  "■  2 

2  e  l+2p-1(2b  +  l) (l-(2b  +  l)  )  (mod  2P ) 


e  1  (mod  2P)  . 


If  m  is  even,  let  m  =  2b  .  Then  (2.5*6)  becomes 


>p-2 


2  e  1+2P  1(2b) (l-(2b) )  (mod  2P) 


e  1  (mod  2P ) . 
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pP-2 

Thus  a  e1  (mod  2P)  for  all  a  . 

For  n  =  2P_^  ,  (2.5*5)  becomes 


2P“3 

a 


1+2P-2m+2P'2(.25~2-l.)m2  +  2P-'3(2P~3-U(-2P-~3.-22(2m)3 


2! 


3! 


! 2P~ 3 ( 2P  3-1)(2p~3-2)(2p~3-^)( 2m) 4 

I  • 

+  higher  powers  of  2  than  2P 
e  l+2p”^m-2p_^m^-2p  ^m^  (mod  2P )  .  (2.5*7) 


If  m  is  even,  i.e.,  if  m  =  2b  ,  (2.5*7)  becomes 


a 


l+2p-1b  (mod  2P )  . 


Then  if  b  is  odd, 


a 


2 


P-3 


1+2P  1  (mod  2P ) 
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or,  if  b  is  even. 


2P-3 

a 


1  (mod  2P)  . 


If  m  is  odd,  i.e.  ,  m  =  2b  +  l  ,  (  2  0  5  ° 7 )  becomes 

a2?  ^  =  l+2p'2(2b+l)-2p-2(2b+l)2-2p“1(2b+l)4  (mod  2P ) 
e  i+2P_1b+2P"2-2P"2-2P_1  (mod  2P ) 

=  l+2p”2(b-l)  (mod  2P)  . 


Then  if  b  is  odd, 


2P-3 


1+2P~’1  (  2c  +  l-l )  (mod  2P  ) 


1  (mod  2P ) 
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or,  if  b  is  even. 


~p-3 

a  e  1+2P  1 ( 2c-l )  (mod  2P ) 

e  1+2P"1  (mod  2P )  . 

From  the  preceding,  it  can  be  readily  seen  that 

a  e  ±3  (mod  8)  . 

Examples  of  multipliers  and  their  resulting  sequences  will 
be  shown  in  the  next  chapter. 

A  variation  to  Lehmer’s  method  consists  of  adding 
a  constant  c  to  equation  (2.5.3)  to  give 

x. ,,  e  a  x.  +  c  (mod  M)  (2.5.8) 

l  +  l  i 

This  variation  is  referred  to  as  the  mixed-congruential 
method.  In  order  to  ensure  the  maximum  period,  conditions 
upon  a,  x,  and  c  may  be  stated  in  the  following  theorem: 


■ 
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Theorem  2.5.2:  The  sequence  defined  by  (2.5.8)  has  full 
period  M  ,  provided  that 


i) 

c 

is 

relatively  prime  to 

M  , 

ii) 

a 

=  1 

(mod  p) 

if  p  is  a 

prime  factor  of 

iii ) 

a 

E  1 

(mod  4) 

if  4  is  a 

factor  of  M  . 

The  proof  of  Theorem  2.5.1  may  be  found  in  Hull  and 
Dobell  (1962). 

The  mi xed-congruential  generators  have  a  few 
advantages  and  some  disadvantages  when  compared  to  the 
multiplicative-congruential  generators.  The  main 
advantage  of  the  mixed  generators  lies  in  the  fact  that 
a  full  period  M  can  be  attained  whereas  it  cannot  be 
attained  with  multiplicative  generators.  Common 
multipliers  for  the  mixed  methods  were  a  =  2P+1  for 
binary  machines  or  a  =  10p+l  for  decimal  machines. 

The  value  of  p  is  the  number  of  binary  bit  positions  or 
the  number  of  decimal  digit  positions  available  to 
represent  an  integer  in  the  computer  being  used.  This 
enabled  a  simple  shift  and  add  rather  than  a  full 
multiplication  to  take  place.  However,  it  was  multipliers 
of  this  type  which  gave  rise  to  rather  poor  statistical 
behaviour  (see  Hull  and  Dobell  (1964)).  With  multiplica¬ 
tive-congruential  generators,  the  statistical  behaviour 
of  sequences  satisfying  the  conditions  of  Theorem  2.5.1 


(q  bom 

: 
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is  generally  acceptable. 

Another  type  of  pseudo-random  number  generator 
which  has  been  proposed  is  an  additive  type.  The 
principal  reason  for  introducing  such  a  generator  is  its 
speed  over  the  multiplicative  methods  since  addition  is 
performed  more  quickly  than  multiplication  in  most 
computers.  One  additive  generator  which  has  been  suggested 
(see  Taussky  and  Todd  (1956))  is  one  involving  a  reduced 
Fibonacci  sequence.  Initially,  F^  =  0  and  F^  =  1  , 
succeeding  numbers  are  generated  by 


F 


n+2 


F  , . +F  (mod  M) 
n+1  n  . 


(2.5.9) 


The  results  of  this  generator  showed  that  the  speed  of 
generation  and  the  period  of  these  numbers  were 
satisfactory,  but  that  the  successive  members  of  F^  were 
not  independent  (see  Taussky  and  Todd  (1956)).  As  a 
remedy,  alternate  members  of  the  sequence  were  chosen. 

The  results  showed  a  more  satisfactory  behaviour,  but 
were  still  not  as  good  as  with  the  multiplicative 
generators.  Green,  Smith  and  Klem  (1959)  suggested  a 
different  method  for  an  additive  pseudo-random  number 
generator . 
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Their  method  used  the  formulae 

X.  =  (X.  ,  +X  .  )  mod  1 

J  J-l  J-n' 

for  decimal  machines,  and 

X.  =  (X.  ,+X.  )  mod  2P 

J  J-l  J -n 

for  binary  machines.  In  this  method,  the  most  recent  n 
random  numbers  must  be  stored.  The  cycle  will,  of  course, 
be  periodic  as  soon  as  the  original  n  numbers  are 
generated.  Tests  performed  by  Green,  Smith  and  Klem 
showed  the  generator  to  be  unsatisfactory.  To  remedy 
this  situation,  alternate  members  were  chosen.  The  results 
from  the  multip licati ve-congruential  generators,  however, 
were  still  more  acceptable. 


CHAPTER  III 


RANDOM  NUMBER  TESTING 


3 . 1  Introduction 

The  purpose  of  this  chapter  is  to  examine  various 
techniques  for  testing  a  sequence  of  digits  or  a  sequence 
of  numbers  for  random  properties.  First,  four  basic 
tests  due  to  Kendall  and  Babingt on-Smith  will  be  examined, 
then  more  recent  tests  specific  to  computer-generated 
pseudo-random  sequences  will  be  discussed.  Some  results 
from  the  tests  are  given  in  the  Appendix  for  the 
multiplicative-congruential  method . 

We  should  note  that  the  sequences  of  pseudo-random 
numbers  were  generated  and  tested  on  an  IBM  7040  computer, 
a  thirty-six  bit  binary  machine.  For  this  reason,  most 
of  the  chapter  will  be  restricted  to  the  consideration  of 
sequences  of  octal  digits.  It  is  a  relatively  simple  task 
to  convert  the  test  to  consider  the  decimal  case. 

3 • 2  Four  Basic  Tests  for  Random  Digits 

In  addition  to  constructing  a  mechanical  random 
digit  generator  (see  2.2),  Kendall  and  Babingt on-Smith 
(1938,1939a)  also  derived  four  basic  tests  for  randomness 
in  a  sequence  of  digits.  The  tests  were  the  frequency 
test,  the  serial  test,  the  poker  test,  and  the  gap  test. 
Each  of  these  tests  will  be  discussed  in  turn. 
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The  frequency  test  is  probably  the  easiest  to 
understand.  In  a  sequence  of  random  octal  digits,  we 
would  expect  to  find  every  digit  occurring  approximately 
an  equal  number  of  times.  For  example,  we  would 
expect  the  octal  digit  5  to  occur  12.5  percent  of  the 
time  in  a  set  of  random  octal  digits.  Any  marked 
departure  from  equality  of  frequencies  would  lead  one 
to  suspect  a  bias  toward  some  digits  and  away  from  others. 
Kendall  and  Babingt on-Smith  used  a  chi-square  test  to 
test  the  hypothesis  of  uniformity,  i.e.  ,  equal  frequencies 
for  each  digit . 

The  serial  test  involves  a  similar  type  of  test 
with  pairs  of  digits.  We  would  expect  that  no  single 
digit  should  tend  to  precede  or  follow  another  digit  if 
a  series  of  digits  were  to  be  locally  random.  In  order 
to  use  the  test,  a  two-way  table  would  be  constructed. 

The  entries  in  the  table  would  be  dependent  on  the 
sequence  of  digits  examined.  For  example,  if  an  octal 
digit  4  followed  an  octal  digit  6  in  the  sequence,  then 
we  would  place  an  entry  in  the  fourth  row  and  sixth 
column  of  the  table.  When  the  entire  sequence  of  digits 
has  been  examined,  the  number  of  entries  in  each  element 
of  the  table  would  be  counted  and  we  would  expect  the 
totals  to  be  approximately  equal.  This  hypothesis  could 
be  tested  using  a  chi-square  test  with  sixty-three  degrees 
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of  freedom. 

As  an  example  of  the  two-way  table  for  serial 
frequencies,  consider  the  sequence  of  digits 

151013324625660211152  . 

Then,  Table  3.2.1  shows  the  serial  frequencies  of  the 
digits . 


0 

1 

2 

3 

4 

5 

6 

7 


0 

0 

1 

0 

0 

0 

0 

1 

0 


1 

1 

2 

2 

0 

0 

1 

0 

0 


2 

1 

0 

0 

1 

0 

1 

1 

0 


3 

0 

1 

0 

1 

0 

0 

0 

0 


456 
0  0  0 
0  2  0 
110 
0  0  0 
0  0  1 
0  0  1 
0  0  1 
0  0  0 


Table  3*2.1  -  Serial  Frequency  of  Some 

Octal  Digits 


7 

0 

0 

0 

0 

0 

0 

0 

0 


A  generalization  of  the  serial  test,  not  suggested 
by  Kendall  and  Babington-Smith ,  would  be  to  compare  two 
digits  which  are  separated  from  each  other  by  n  digits, 
where  n  is  any  positive  integer.  We  would  expect  the 
results  from  this  modification  to  be  exactly  the  same 
for  any  value  of  n  . 
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The  poker  test  examines  groups  of  five  digits  and 
tests  their  value  as  a  poker  hand.  The  events  and  their 
expected  occurrance  are  shown  in  Table  3.2.2. 


Event 

Description 

Example 

Probability 

Bust 

All  digits  different 

abode 

.205 

1  Pair 

Two  digits  the  same 

aabcd 

.513 

2  Pairs 

Two  pairs  of  digits 

aabbc 

.153 

the  same 

Triple 

Three  digits  the  same 

aaabc 

.103 

Full  House 

One  pair,  one  triple 

aaabb 

.017 

4  of  a  Kind 

Four  digits  the  same 

aaaab 

.009 

5  of  a  Kind 

All  digits  the  same 

aaaaa 

.  000 

1.000 

Table 

3.2.2  -  Probabilities 

for  Selected 

Poker  Hands 

In  order  to  use  the  poker  test,  the  sequence  of  digits 
is  divided  into  groups  of  five  digits  and  the  observed 
frequency  of  the  possible  poker  hands  is  compared  to  the 
theoretical  frequency  by  a  chi-square  test.  Since  the 
probabilities  of  a  full  house,  four  of  a  kind  and  five 
of  a  kind  are  relatively  small,  the  frequencies  of  these 
events  are  usually  grouped  together. 

The  gap  test  compares  the  expectations  with  the 
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frequency  of  the  gaps  between  two  successive  equal  digits. 
A  gap  may  be  defined  as  follows.  If  a  digit  n  is 
separated  from  the  same  digit  n  by  m  digits,  this 
constitutes  a  gap  of  length  m  for  the  digit  n  .  The 
probabilities  for  octal  digit  gaps  may  be  summarized 
in  Table  3.2.3. 
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Length  of  Gap 
0 
1 
2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16-2-0 

21-25 

>26 

Table  3.2.3  -  Gap 


Probability 
.  125 
.109 
.095 
.083 
.073 
.064 
.056 
.049 
.043 
.038 
.033 
.029 
.025 
.022 
.019 
.017 
.058 
.030 
.032 

1.000 

Probabilities  for 


Octal  Digits 


>- 
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For  the  gap  test,  we  may  test  all  eight  octal  digits. 
However,  usually  it  is  sufficient  to  test,  for  example, 
the  gaps  for  the  octal  digit  2  since  we  would  expect 
similar  results  for  any  other  octal  digit. 

Kendall  and  Babington-Smith  were  careful  to 
caution  that  if  two  locally  random  sets  of  digits, 
i.e.,  sets  whose  digits  pass  the  preceding  four  tests, 
were  combined,  the  resulting  sequence  of  digits  would 
not  necessarily  be  locally  random.  Furthermore,  as  the 
number  of  random  digits  in  the  set  increased,  the  number 
of  bad  patches  appearing  without  local  randomness  was 
bound  to  increase. 

The  chi-square  test  was  used  by  Kendall  and 
Babington-Smith  to  measure  permissable  deviation  from 
expectations.  The  table  of  random  digits  produced  by 
the  two  authors  satisfied  the  four  tests  and  was  widely 
used  until  the  production  of  RAND’s  million  random  digits. 

3 . 3  Modifications  to  Kendall  and  Babington-Smi th ’ s  Tests 

The  frequency  test  can  be  applied  in  two  different 
ways.  First,  a  count  can  be  made  of  the  frequency  of 
occurrence  of  each  octal  digit  in  the  set  of  pseudo-random 
numbers.  Second,  a  count  can  be  made  of  each  digit 
position  of  the  pseudo-random  number.  For  example,  let 
us  consider  the  ten  sequences  of  octal  digits  in  Table 


3.3*1  which  were  generated  by  the  relation 

xi  +  l  E  (21^  +  3)(mod  2^)  ,  where  xQ  =  (123456789) 

366017263477 

073715432675 

101107120467 

353560761645 

274045325357 

236747600315 

235035601147 

227614603465 

211100612637 

140621640335 

Table  3.3*1  -  First  Ten  Numbers  from 

xi+1  =  ( 2l8+3) x1 (mod  2^) 

x  0  =  ( 123456789 ) 1Q 


The  first  frequency  table  is  shown  in  Table  3.3*2 
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Digit 

0 

1 

2 

3 

4 

5 

6 

7 


Frequency 
1 6 

17 
12 
1 6 
12 
13 

18 
16 

120 


Table  3.3.2  -  Observed  Frequencies  of  Octal 

Digits  of  Table  3.3.1 


The  frequency  of  the  octal  digits  would  then  be  compared 
to  the  expected  frequencies  of  fifteen  by  a  chi-square 
test  with  seven  degrees  of  freedom. 

The  second  frequency  test  mentioned  above  would 
be  tallied  in  a  similar  fashion  except  that  there  would 
be  twelve  separate  tables  each  with  eight  entries.  An 
example  of  this  type  of  table  is  illustrated  in  Table 


3.3.3 


36 


igit 

0 

1 

2 

3 

Octal  Position 

4  5  6  7 

8 

9 

10 

11 

0 

1 

1 

1 

3 

2 

2 

0 

3 

3 

0 

0 

0 

1 

2 

1 

2 

2 

3 

1 

1 

1 

2 

1 

1 

0 

2 

5 

1 

0 

0 

1 

0 

1 

2 

2 

0 

0 

0 

3 

2 

2 

2 

0 

1 

0 

1 

1 

2 

3 

2 

0 

H 

0 

1 

1 

0 

2 

1 

1 

1 

0 

3 

2 

0 

5 

0 

1 

1 

1 

0 

3 

0 

0 

1 

0 

1 

5 

6 

0 

1 

2 

2 

1 

0 

5 

2 

0 

3 

2 

0 

7 

0 

2 

1 

2 

0 

3 

1 

0 

0 

0 

2 

5 

Table  3° 3- 3  -  Observed  Frequencies  for  Individual 

Octal  Digits  from  Table  3.3*1 

Of  course,  one  can  see  immediately  from  inspection 
of  Table  3.3*1  that  the  last  digit  will  always  be  odd0 
This  is  a  result  of  the  restrictions  on  and  a  in 

the  multiplicative  congruential  method,  The  first  digits 
appearing  in  Table  3*3.1  are  always  0,  1,  2,  or  3  because 
of  the  octal  representation  of  numbers  in  the  7040  -  the 
first  octal  digit  always  contains  the  sign  of  the  numbers. 
These  phenomena  force  us  to  neglect  some  of  the  octal 
digit  positions  in  a  pseudo-random  number. 

A  further  variation  to  the  second  method  is  to 
divide  the  interval  (0,1)  into  equal  subintervals  and  to 
test  the  observed  frequency  of  numbers  within  the  interval 
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with  the  expected  frequency .  However,  the  number  of 
intervals,  usually  restricts  the  test  to  the  first  few 
octal  digits .  For  example,  if  we  divide  the  interval 
into  eight  equal  segments,  then  we  would  only  be  testing 
the  first  octal  digit. 

The  serial  test  for  pseudo-random  numbers  can  be 
performed  in  a  manner  similar  to  that  of  the  frequency 
test.  Let  us  consider  the  sequence  of  numbers  generated 
in  Table  3.3.1.  Then  the  serial  test  performed  on  the 
entire  set  of  numbers  would  give  the  results  shown  in 
Table  3.3°40 


0 

1 

2 

3 

4 

5 

6 

7 


0 


1 


7 


2 

3 

1 

0 

3 

1 

5 

1 


3 

4 
2 
1 
0 
1 
3 
3 


0 

2 

1 

2 

0 

3 

1 

3 


4 

0 

2 

1 

1 

3 

2 

2 


2 

3 

0 

2 

0 

1 

2 

2 


0 

2 

1 

6 

2 

0 

1 

1 


2 

2 

3 

2 

3 

2 

1 

3 


3 

1 

2 

2 

3 

1 

3 

1 


Table  3 ° 3  -  4  -  Serial  Frequencies  for  the  Octal 

Digits  in  Table  3*3.1 


The  observed  frequencies  are  compared  to  the  expected 
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frequencies  by  means  of  the  chi-square  test. 

A  second  method  for  calculating  serial  frequency 

would  be  to  examine  each  octal  digit  position  separately 

in  a  manner  similar  to  that  of  the  frequency  test. 

Hull  and  Dobell  (1964)  referred  to  a  paper  by  Good 

(1953)  in  which  the  results  from  the  frequency  test  and 

the  serial  test  were  combined  to  test  a  sequence  for 

randomness.  It  was  shown  in  the  latter  paper  that  if 
2 

represents  the  chi-square  from  the  frequency  test 
2 

and  X2  represents  the  chi-square  from  the  serial  test, 

2  2  2 

then  =  x2  ”  X1  s  asymP^0^ca^y  distributed  as  a 

2 

chi-square  distribution  with  v  -  v  degrees  of  freedom 
where  v  =  8  for  the  case  of  octal  numbers.  If  we  are 
examining  the  frequency  of  numbers  within  equal 
subintervals,  then  v  is  the  number  of  intervals. 

The  poker  test  can  also  be  formulated  in  two  ways. 
First,  we  could  isolate  the  digits  of  a  pseudo-random 
number  serially,  five  at  a  time,  and  test  the  observed 
frequency  to  the  expected  frequency.  Second,  we  could 
examine  the  individual  octal  positions  in  the  same 
manner  as  in  the  frequency  test  and  perform  a  similar 
chi-square  test. 

The  gap  test  could  also  be  treated  in  the  same 
manner,  but  its  use  seems  to  have  been  emphasized  less  in  the 


* 
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current  literature  than  the  use  of  the  serial  and 
frequency  tests. 

A  further  application  of  the  chi-square  test  from 
the  values  of  chi-square  obtained  from  Kendall  and 
Babingt on-Smith  tests  was  suggested  by  Green,  Smith,  and 
Klem  (1959).  When  the  values  of  chi-square  are  calculated 
for  a  sufficiently  large  number  of  frequency  and  serial 
tests,  we  should  obtain  a  whole  range  of  chi-square 
values o  The  frequency  of  these  values  should  in  turn, 
also,  satisfy  a  chi-square  distribution  with  n  intervals 
between  0  and  100  percent.  This  additional  test  has  not 
been  used  to  a  wide  extent  since  it  requires  the  generation 
of  several  blocks  of  pseudo-random  numbers  and  the  testing 
of  each  individual  block. 

3  0  4  Other  Tests  for  Randomness 

Although  Kendall  and  Babingt on-Smith ' s  tests  still 
form  the  basis  for  testing  pseudo-random  numbers,  other 
tests  have  been  suggested  in  recent  years. 

Gruenberger  and  Mark  (1951)  suggested  a  test  for 
use  in  Monte  Carlo  calculations.  Since  many  Monte  Carlo 
techniques  require  the  use  of  special  coordinates,  two 
successive  pseudo-random  numbers  can  be  used  to  determine 
( x ,y ) -coordinates  within  the  unit  square.  The  test 
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suggested  an  examination  of  the  square  of  the  distance 

between  two  successive  points.  For  two  points,  the 

probability  that  the  square  of  the  distance  between  them 
2 

is  a  is  given  by 


P 


2 

tt  a 


8a3 

3 


+ 


for  0 


2 

<  a 


<  1.0 


and 


P  =  —  + 

3 


77 


-2)a2+4(a2-l)1/2  +  |(a2-l) 


3/2 


4 


2  -1 
a  sec  a 


for  1.0 


2 

<  a 


<  2.0 


(3.4.1) 


In  order  to  perform  the  test,  the  square  of  the  distance 

is  computed  and  the  distribution  of  the  results  is 

compared  to  the  theoretical  by  the  chi-square  test.  The 

2 

cumulative  probabilities  for  discreet  a  are  shown  in 
Table  3.4.1  (see  Gruenberger  and  Mark  (1951))  as  well  as 
the  probabilities  for  each  of  the  intervals  0.0  (0.1) 


2.0. 


4l 


a 

P ( ( a-0 . 1 )  <d 

0.1 

.234832 

0.2 

.174973 

0.3 

.139495 

0.4 

. 112718 

0.5 

.090971 

0 . 6 

.072614 

0.7 

.056748 

0.8 

.042814 

0.9 

.030430 

1.0 

.019333 

1.1 

.010777 

1.2 

.006345 

1.3 

.003740 

1.4 

.002138 

1.5 

.001154 

1.6 

.000572 

1.7 

.000246 

1.8 

.000084 

1.9 

.000017 

2.0 

.000001 

Table 

3.4.1  -  Probabilit; 

Cumulative  Probability 
.234832 
.  409805 
.549300 
.662018 
.752987 
.  825601 
.882349 
.925163 
.955593 
.974926 
.985703 
.992048 
.995788 
.997926 
.999080 
.999652 
.999898 
.999982 
.999999 
1.000000 

Associated 


with  the  d  -Test 


. 

Gorenstein  (1967)  described  tests  on  the  moments 
of  the  uniform  distribution.  He  computed  the  mean  and 
the  second  and  third  moments  for  equal  blocks  of 
generated  pseudo-random  sequences.  The  expected  values 
are  shown  in  Table  3*4.2  for  a  uniform  generator  with 
modulus  2P . 

Moment  Distribution  Expected  Value 


First 


1 

2PN 


r  N 
1  i  =  l 


1/2 


Second 


1 

2P2PN 


I 


N 

i  =  l 


1/3 


Third 


_ 1 _  y  n  x3 

2P2P2PN  1=1  1 


1/4 


Table  3*4.2  -  Moments  for  Gorenstein’s  Tests 


McLaren  and  Marsaglia  (1965) 
applications  of  order  statistics  to 
tests  for  pseudo-random  sequences, 
based  on  the  following  properties. 


used  two  simple 
perform  additional 
Their  tests  were 
From  a  set  of  n 


pseudo-random  numbers,  a  maximum  element,  Max,  and  a 


iiamoni 


3- 


4  3 


minimum  element,  Min,  are  chosen.  If  the  n  pseudo-random 
numbers  {x^,...,x  }  are  uniformly  distributed,  then  Max 
should  have  a  distribution  P(Max  <  a)  =  F(a)  =  an  for 
0  <  a  <  1  .  Similarly,  Min  should  have  a  distribution 
P(Min  >  a)  =  $(a)  =  (l-a)n  .  The  chi-square  test  is 
again  used  to  test  equal  subintervals  of  a  large  block  of 
pseudo-random  numbers. 

3 . 5  Additional  Comments 

Peach  (1961)  discovered  an  interesting  phenomenon 
about  the  mixed  congruential  methods.  He  considered  the 
example  xn+]_  =  9xn  +  13  (mod  32)  which  generated  the 
sequence 


0 

13 

2 

31 

4 

17 

6 

3 

8 

21 

10 

7 

12 

25 

14 

11 

16 

29 

18 

15 

20 

1 

22 

19 

24 

5 

26 

23 

28 

9 

30 

27 

Each  of  the  above  rows  comprises  a  quarter-period  of  the 
sequence.  Upon  close  examination  of  the  half-periods, 
it  can  be  seen  that  each  corresponding  number  differs  from 
the  other  by  exactly  sixteen.  Similarly,  for  each 
quarter-period,  the  difference  is  eight,  and  for  each 
eighth-period  the  difference  is  four.  In  general,  it  was 
observed  that  half-periods  differed  by  2  ”  ,  quarter- 


>- 
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periods  differed  by  2m  ^  ,  eighth-periods  differed  by 
rn  q 

2  ,  etc.  These  periodicities  lead  us  to  caution 

that  pseudo-random  numbers  generated  by  the  mixed- 
congruential  method  contain  patterns  and  periodicities 
which  act  as  constraints  upon  their  variability. 

With  the  preceding  example  in  mind,  let  us  now 
examine  a  multiplicat i ve-congruential  generator, 


cally 

xi  +  l 

29x. 

1 

(mod 

64) 

with  x0  =  11 

uence 

generate 

?d  is 

11 

63  35 

55 

59 

47 

19  39 

43 

31  3 

23 

27 

15 

51  7  11 

Here,  we  have  the  same  phenomenon  as  before  -  half-sequences 
differ  by  thirty-two,  quarter-sequences  differ  by  sixteen. 
Greenberger  (1965)  discovered  that  the 

l8  88 

multiplicative  congruential  generator  x1  +  1  E  (  210  +  3 )  (mod  2^')) 

produced  a  second-order  serial  correlation.  In  order  to 

rectify  the  situation,  a  series  of  five  successive 

1 8 

multipliers  close  to,  but  not  equal  to,  2±  +3  were 

used,  the  successive  multipliers  being  chosen  "at  random". 

In  the  multiplicative-congruential  method,  only 
half  of  the  odd  numbers  can  possibly  occur.  For  example, 
if  a  e  3  (mod  8)  and  Xq  e  1  (mod  8)  ,  then 
x^  +  ^  =  1  (mod  8)  or  xj_  +  i  E  3  (mod  8)  .  In  order  to  obtain 
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the  other  odd  numbers,  we  would  require  a  further 
multiplier  of  5  (mod  8).  Greenberger  used  five 
multipliers  which  should  scramble  the  sequence  even 
further . 

Hutchinson  (1966)  suggested  a  generator  which 
revived  Lehmer’s  original  proposal.  For  a  machine 
capable  of  representing  an  integer  of  maximum  absolute 
value  2P  -  1  ,  the  largest  prime  less  than  2P  is 
used  as  the  modulus.  The  value  of  a  chosen  is  a 
primitive  root  of  the  prime.  For  the  7040,  the  modulus 
is  2 ^  -  31  j  a  is  5^  or  5^  .  For  the  360,  the 
modulus  is  2^  -  1  ,  a  is  2^  +  1  .  The  distinct 
advantage  of  the  Hutchinson  method  when  compared  to 
the  more  commonly  used  mult iplicat i ve-congruent ial 
method  where  M  =  2^  is  that  the  least  significant 
bits  appear  to  be  as  random  as  the  most  significant 
bits.  In  order  to  understand  the  generator,  let  us  look 
at  the  sequence  generated  by  a  =  24  ,  Xq  =  11  ,  modulus 
31  .  We  can  assume  that  the  word  size  for  this  example 

5 

could  accommodate  a  number  as  large  as  2  - 

sequence  generated  would  be 


1  .  The 


'  - 
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11 

16 

12 

9 

30 

7 

13 

2 

17 

5 

27 

28 

21 

8 

6 

20 

15 

19 

22 

1 

24 

18 

29 

14 

26 

4 

3 

10 

23 

25 

11 

If  we  examine  the  differences  between  successive  elements 
of  the  series,  we  notice  that  the  differences  are 

5  27  28  21  8  6 

This  sequence  is  exactly  the  same  sequence  of  numbers 
generated  except  that  the  starting  value  has  changed. 
Similarly,  if  we  examine  the  differences  of  lag  2,  the 
same  phenomenon  occurs.  In  Table  3.4.3,  the  lag  is 
specified  as  well  as  the  first  few  terms  of  the  sequence 


of  differences. 


>• 
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Lag  Sequence 


1 

5 

27 

28 

21 

2 

1 

24 

18 

29 

3 

29 

14 

26 

4 

4 

19 

22 

1 

24 

5 

27 

28 

21 

8 

6 

2 

17 

5 

27 

7 

22 

1 

24 

18 

8 

6 

20 

15 

19 

9 

25 

11 

16 

12 

10 

16 

12 

9 

30 

11 

17 

5 

27 

28 

12 

10 

23 

25 

11 

13 

28 

21 

8 

6 

14 

26 

4 

3 

10 

15 

9 

30 

7 

13 

16 

4 

3 

10 

23 

17 

8 

6 

20 

15 

18 

11 

16 

12 

9 

19 

21 

8 

6 

20 

20 

13 

2 

17 

5 

21 

7 

13 

2 

17 

22 

18 

29 

14 

26 

23 

3 

10 

23 

25 

24 

15 

19 

22 

1 

25 

24 

18 

29 

14 

26 

23 

25 

11 

16 

27 

30 

7 

13 

2 

28 

12 

9 

30 

7 

29 

14 

26 

4 

3 

Table  3.4.3 

-  First  Four 

Members 

of  Sequences 

with  Lags  up  to  29 
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The  significance  of  the  results  of  Table  3.4.3 
are  not  yet  clear  and  might  prove  interesting  for  further 
study . 

3 • 6  Concluding  Remarks 

It  would  be  worth  stating  some  pertinent  remarks 
concerning  pseudo-random  number  generators.  Hull  and 
Dobell  (1964)  point  out  that  any  generators  which  pass 
certain  tests  are  not  completely  guaranteed  of  being 
acceptable.  Some  sequences  of  numbers  may  be  perfectly 
acceptable  as  far  as  the  tests  are  concerned,  but  will 
fail  in  the  particular  application  for  which  they  are 
being  used.  , It  is  imperative  that  we  be  able  to  at  least 
have  a  general  idea  of  the  expected  answer  before 
attempting  the  use  of  pseudo-random  numbers.  Lehmer's 
definition  of  a  pseudo-random  sequence  as  given  by 
Taussky  and  Todd  (1955)  is:  "A  vague  notion  embodying 
the  idea  of  a  sequence  in  which  each  term  is  unpredictable 
to  the  initiated  and  whose  digits  pass  a  certain  number  of 
tests,  traditional  with  statisticians  and  depending 
somewhat  on  the  uses  to  which  the  sequence  is  to  be  put". 


CHAPTER  IV 


THE  GENERATION  OF  NON-UNIFORMLY 
DISTRIBUTED  RANDOM  VARIABLES 

4 . 1  Introduction 

The  main  purpose  of  this  chapter  is  to  examine 
methods  for  converting  a  uniformly  distributed  random 
variable  to  a  normally  distributed  random  variable.  In 
addition,  two  other  distributions  of  interest  -  the 
exponential  distribution  and  the  Poisson  distribution 
will  be  considered. 

First  of  all,  however,  let  us  consider  the 
transformation  of  a  given  random  variable  £  with  a 
known  probability  density  function  to  a  random  variable 
n  with  another  probability  density  function.  The 
following  well  known  theorem  from  Freund  (1962)  is 
fundamental  to  the  material  discussed  in  this  chapter. 

Theorem  4,1,1:  If  the  probability  density  function  of 
xeC  is  given  by  f(x)  and  the  function  given  by 
y  =  h(x)  is  differentiable  and  either  increasing  or 
decreasing  for  all  values  within  the  range  of  £  ,  then 
the  probability  density  function  of  yen  is  given  by 


g(y )  =  f(x) 


dx 

dy 


j 


(4.1.1) 
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where  -7—  i-  0  . 
dy 

Proof:  Let  F  (a)  represent  the  value  of  the  distribution 
function  of  £  at  a  .  Then,  the  probability  that  5 
assumes  a  value  less  than  or  equal  to  a  is 

a 

F  (a)  =  P(£<a)  =  /  f(x)dx  (4.1.2) 

—  00 

We  shall  assume  that  y  =  h(x)  is  increasing.  (The 
proof  may  be  modified  when  h(x)  is  decreasing).  Then, 

F  (a)  also  represents  the  probability  that  n  assumes 
a  value  less  than  or  equal  to  h(a)  ,  i.e., 

a 

F  (h(a))  =  /  f(x)dx  .  (4.1.3) 

^  Loo 


If  the  change  of  variable  y  =  h(x)  ,  x  =  H(y)  is 
performed  in  (4.1.3),  then 


h  ( a) 

P  (h(a) )  =  /  f(H(y ) )  H'(y)dy  ( ^ . 1 . 4 ) 

n 

_  rv \ 
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for  any  real  number  h(a)  within  the  range  of  n  . 
Therefore,  if  H’(y)  exists,  the  integrand  in  equation 
(4.1.4)  gives  the  probability  density  function  of  n  . 
That  is,  the  probability  density  function  of  n  is 
given  by 


g(y )  =  f (H(y ) )  H» (y )  , 


1 

sHTyT 


¥  0 


(4.1.5) 


which  can  be  written 


g(y)  -  f(*>  § 

"1  r\  y 

since  =  0  ,  H(y)  =  x  and  HT(y)  =  ^  . 

Keeping  (1962)  shows  that  if  F(x)  is  the  distribution 
function  of  X  ,  and  if  the  transformation  Y  =  F(X)  is 
made,  then  Y  has  a  uniform  distribution  on  (0,1).  In 
4.3,  we  shall  see  how  the  preceding  theory  can  be  used 
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to  derive  the  exponential  distribution. 

4 c  2  The  Normal  Distribution 

Since  normally  distributed  random  variables  are 
important  in  statistics,  conversions  from  a  uniformly 
distributed  random  variable  to  a  normally  distributed 
random  variable  have  been  studied  and  reviewed  extensively 
in  the  literature.  Among  the  most  important  conversion 
techniques  are  those  considered  by  Box  and  Muller  (1958), 
Hastings  (1955),  Juncosa  (1953),  Marsaglia  (1964),  and 
von  Neumann  (1951). 

Von  Neumann’s  technique  depends  on  an  acceptance- 
rejection  procedure.  In  order  to  generate  normal 
deviates  X  in  the  region  -b  <_  X  <_  b  ,  two  uniform 
random  deviates  and  are  generated.  Then 

Y  =  -2b2(U1~.5)^  is  computed.  If  logg  £  Y  ,  then 
X  =  b(2U1~l)  is  used  as  a  normal  deviate.  Otherwise, 
if  logg  u2  >  Y  ,  then  the  pair  (U^U^)  is  rejected  and 
another  pair  of  uniform  random  deviates  would  be 
generated  and  the  procedure  described  above  is  repeated. 

The  von  Neumann  approach  is  inefficient  for 
generating  normal  deviates  especially  when  values  beyond 
three  standard  deviations  from  the  mean  are  needed.  The 
inefficiency  is  recognized  as  soon  as  an  expression  is 
found  for  the  probability  that  a  pair  (U^jU^)  will  be 


used  to  generate  a  normal  deviate.  The  probability  may 
be  expressed  as 
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p  { U2  1 


e-2b2(U1-.5)2} 


=  7 


1 


-2b<:(U1-.5) 


dU. 


0 


(4.2.1) 


which  asymptotically  approaches  (1/b)  A/2  »  If  b 

becomes  small,  then  the  probability  increases  proportionately. 

The  Hasting’s  approach  uses  a  Pade  approximation 
to  transform  a  uniform  random  deviate  to  a  normal  random 
deviate.  The  most  accurate  method  proposed  obtains  a 
normal  random  deviate  X  from  a  uniform  random  deviate 
q  using 


an+an  n+a 9n 

X  =  X#(q)  =  n  -  ( — - - - p)  (4.2.2) 

l+bp+b^n+bpri 


where 


aQ  =  2.515517 
a;L  =  0.802853 
a2  =  0.010328 


bQ  =  1.432788 
b1  =  0.189269 
b2  =  0.001308 


>§3|< 
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n 


(l/q)2 


and 


q 


/ 

xtq) 


dt 


0  _<  q  <_  .  5 

If  q  lies  in  the  interval,  .5  <  q  <_  1  ,  then  a 
transformation  must  be  performed  on  q  to  reduce  it 
to  the  [0,5].  The  transformation  used  is 


q  =  1  -  q  e 

The  value  of  X  is  found  from  (4.2.2)  and  the  normal 
random  deviate  used  is  the  value  of  (4.2.2)  with  the  sign 
changed . 

Hasting’s  approach  gives  results  which  are  accurate 
to  within  +  4  x  10”^  of  the  correct  result.  However, 
the  time  taken  to  generate  normal  deviates  is  greater 
than  for  some  other  methods.  In  order  to  compensate  for 


>- 
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this,  Hasting's  derived  another  Pade  approximation  using 
fewer  constants .  The  conversion  is  given  by 


X*(q)  =  n 


an+a, n 

{ - - — - — 


1+bQn+b-^n 


2 


(4.2.3) 


where  aQ  =  2.30753  bQ  =  0.99229 

a1  =  0.27061  b1  =  0.04481 

and  q  and  n  are  given  by  equation  (4.2.2).  The 
decrease  in  time  is  compensated  by  the  difference  in  the 
maximum  absolute  error.  In  (4.2.3),  the  maximum  absolute 
error  is  3  x  10  . 

Juncosa's  approach  makes  use  of  the  Central  Limit 
Theorem  which  may  be  stated  as  follows  (Keeping  (1962)). 

Theorem  4.2.1:  Let  X^,  ,  X^,...,  X^  be  independent 

random  variates  all  having  the  same  distribution  with 

2 

mean  y  and  variance  a  ,  but  not  necessarily  normal. 

Let  the  standarized  variate  corresponding  to  X.  be 

J 


-  y 


5 


X. 


a 


(4.2.4) 
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and  let  Yn  be  defined  by 


n 

l  z 

V  =  J=1  ' 

n  /n 


=  /n  Z  , 


(4.2.5) 


where  Z  is  the  arithmetic  mean  of  the  Z.  .  Then,  as 

J 

n  ^  00  ,  Y  tends  to  a  standard  normal  variate. 

5  n 

For  the  uniform  distribution  on  (0,1),  it  can  be  s'hown  that 


v  =  1/2 


and 


a  =  1/12 


Therefore,  from  (4.2.4), 


=  (X .- . 5)  /12 
J 


Z. 

3 


(4.2.6) 
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and 


l 

i=i 


(X  -,5)/l2 

•V 


(4.2.7) 


The  Juncosa  approach  is  reasonably  fast  and  requires 
little  memory  space.  However,  for  every  random  normal 
deviate  generated,  a  sufficiently  large  number  of  random 
uniform  deviates  must  be  available.  In  Juncosa’ s  report, 
sixty-four  uniform  random  numbers  were  generated  for 
each  normal  variate.  The  Appendix  shows  plots  for  some 
values  of  n  . 

Marsaglia’s  method  is  particularly  fast.  The  normal 
distribution  is  expressed  as  a  linear  combination  of 
three  separate  functions, 

g ( x )  =  0.9578gl(x)  +  0.0395g2(x)  +  0.0027g3(x)  (4.2.8) 


g(x) 


2 


/2  TT 


where 


3- 

' 
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The  distributions  of  g(x),  g^(x),  g£(x)  and  g  (x) 
shown  in  Figures  4.2.4,  4.2.1,  4.2.2,  and  4.2.3 
respectively.  (See  Marsaglia  (1964)). 


Figure  4.2.1  -  g-^(x) 


Figure  4.2.2  -  g^(x) 


are 


59 


Figure  4.2.3  -  g^(x) 


Figure  4.2.4  -  g( x ) = . 9578g]_ ( x ) + . 0395g2 ( x ) + . 0027g3 ( x ) 
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The  function  g^(x)  is  very  simple  to  evaluate  and  is 
used  ninety-six  to  ninety-seven  percent  of  the  time.  The 
other  three  or  four  percent  of  the  time  a  modification 
involving  a  complicated  set  of  instructions  is  used  to 
make  the  function  agree  with  the  normal  distribution. 

The  flow  chart  for  a  binary  machine  is  shown  in  the 
Appendix. 

The  Marsaglia  method  is  fast  since  it  consists 
mainly  of  a  table  look-up.  However,  it  requires  a  large 
amount  of  storage  to  store  all  the  constants  necessary. 
The  method  could  be  made  to  generate  the  entire  range  of 
values  if  a  plus  or  a  minus  sign  were  placed  randomly 
in  the  result 

The  Box  and  Muller  approach  converts  two  uniform 
random  deviates  to  two  normal  random  deviates  directly. 
The  method  takes  the  two  uniform  random  deviates 
and  from  the  interval  (0,1)  and  transforms  them  as 

follows : 


X-^  =  (-2  logeU-^)^  cos  2-nU^ 


X2  = 


(-2  log  U-  ) ^  sin  2ttU2 


(4.2.9) 
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(X^jX^)  will  then  be  a  pair  of  independent  normal  random 
variables  with  mean  0  and  variance  1. 

In  order  to  justify  the  above,  from  ( 4 . 2 . 9  )  ,  the 
inverse  relationships  are 


U 


1 


U 


2 


1 

arctan 
2  IT 


2 
X 


(4.2.10) 


It  can  be  shown  that  the  joint  density  of  X-^,  X^  is 


f(x1,x2) 


/2  IT  /2  TT 


=  f(X1)f(X2)  . 


(4.2.11) 


which  gives  (X^,X2)  as  two  independent  normally  distributed 
random  deviates  with  mean  0  and  variance  1. 
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The  direct  approach  of  Box  and  Muller  gives  better 
accuracy  than  comparable  methods  and  involves  little 
memory  space.  An  added  advantage  to  this  method  is  that 
values  beyond  three  standard  deviations  of  the  mean  are 
reliable.  Furthermore,  the  transformation  from  uniform 
to  normal  is  exact  and  every  uniform  variable  can  be 
t rans  formed . 


4 . 3  The  Exponential  Distribution 

From  4.1,  a  method  can  be  derived  to  convert  a 
uniform  random  variable  E,  to  an  exponentially  distributed 
random  variable  n  „ 

Let  f(x)  be  the  probability  density  function  of 
the  uniformly  distributed  random  variable  £  ,  and  let 
h(y)  be  the  probability  density  function  of  the 
exponentially  distributed  random  variable  n  .  Then 


h(y )  =  Xe  Xy 


(4.3.1) 


Let  represent  a  uniformly  distributed  pseudo-random 

number  from  which  we  wish  to  generate  ,  an 

exponentially  distributed  random  number.  From  4.1.2  and 
4.1.4 
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W 


/ 


n 


i 


—  oo 


Ae 


-Ay 


dy 


(4.3.2) 


-e 


An 


i  +  1 


There  fore 


1  - 


and 


n±  =  j  loge  (l-5±)  (4.3.3) 

If  the  numbers  ^  are  uniformly  distributed  on  (0,1), 
then,  by  solving  (4.3*3),  it  is  possible  to  compute  a 
sequence  of  random  variables  n^  with  the  exponential 


distribution . 
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4  o  4  Poisson  Distribution 

We  will  now  turn  our  attention  to  the  generation 
of  a  discrete  random  variable  distribution  by  considering 
the  Poisson  distribution.  The  probability  distribution 
for  Poisson  is 


x  -X 

f(x,A)  =  — |t~ 


(4.4.1) 


where 


x  =  0 , 1  j  2  , .  .  . 


and 


A  =  n6  ,  (n-*°°,e+0)  ,  is  constant 


The  cumulative  probability  density  distribution  is  given  by 


F(x,A) 


x 


-  I 

i  =  0 


e 


-A 


i ! 


(4.4.2) 


. 
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In  order  to  generate  random  variables  satisfying  a  Poisson 
distribution  from  a  continuous  uniform  distribution,  we 
would  proceed  in  the  following  manner.  We  would  generate 
a  uniform  deviate  FL  in  (0,1).  If 


( 

where 


k 


I 

i  =  0 


X 


TT 


3 


then  the  value  j  is  used  as  the  Poisson  deviate.  This 
procedure  would  be  repeated  for  a  sufficient  number  of 
trials.  The  distribution  of  the  j's  would  then  satisfy 
the  Poisson  distribution  for  the  particular  X  chosen. 
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CHAPTER  V 


SOME  SIMPLE  APPLICATIONS 

5 • 1  Buff on  * s  Needle 

A  rough  approximation  to  the  value  tt  may  be 
determined  by  considering  Buffon’s  needle  problem  (see 
Keeping  (1962)).  The  problem  consists  of  a  needle  of 
length  a  thrown  onto  a  table  which  has  been  ruled  with 
equidistant  parallel  lines  a  distance  a  ( £<a)  apart. 
Buff on,  the  celebrated  French  naturalist  posed  the 
problem:  What  is  the  probability  that  the  needle  will 

Intersect  one  of  the  lines?  We  may  form  an  empirical 
formula  for  the  probability  as  follows.  If  we  take  the 
x-axis  as  being  one  of  the  parallel  lines,  then  the 
x-coordinate  of  the  centre  of  the  needle  will  be 
immaterial,  i.e. ,  we  need  only  consider  the  y-coordinate 
of  the  centre  of  the  needle  and  the  angle  0  which  the 
needle  makes  with  the  x-axis.  The  needle  will  cross  a 
parallel  straight  line  if  the  distance  y  from  the 
centre  of  the  needle  to  the  nearest  line  is  less  than  or 
equal  to  \  l  sin  0  .  (See  Figure  5.1.1). 


, 
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The  probability  that  the  needle  will  cross  the  nearest 
straight  line  is  given  by 


P 


r tt  rh  l  sin  0  N 

J  j  f (0,y )  dy  d0 

0  0 


r  / 
0  0 


a/2 


f(0,y)  dy  d0 


(5.1.D 


Since  we  have  assumed  that  the  needle  has  been  tossed  at 
random,  then  all  values  of  0  and  y  are  equally  likely 
and,  as  a  result,  f(0,y)  is  constant  and  (5.1.1)  becomes 
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P  = 


("ir  i  sin  0  , 

Ji  dy  d0 

0  0 


jtt  j-a/2 
0  0 


dy  d0 


(5.1.2) 


=  iA 

tt  a 


(5.1.3) 


This  experiment  was  performed  in  the  manner  described 
below.  Two  random  numbers,  y  ,  and  y^  ,  were  generated 
within  the  interval  (0,1).  The  centre  of  the  needle  was 
assigned  to  coordinate  y^  .  In  order  to  calculate  the 
value  of  0  ,  y^  was  converted  to  radians  by  setting 
0  =  tt y 2  .  The  test  to  compare  y^  to  \  i  sin  0  was 
performed.  If  y^  >  \  l  sin  0  ,  then  the  total  number  of 
throws  was  incremented  by  one.  If  y^  <_  \  l  sin  0  ,  then 
the  total  number  of  throws  was  incremented  by  one  and  the 
total  number  of  times  a  line  was  crossed  was  incremented 
by  one.  When  a  sufficiently  large  number  of  throws  was 
performed,  the  ratio,  p  ,  of  crosses  to  the  total  tosses 
was  computed.  The  value  of  u  was  approximated  using 


TT 


2% 

pa 


(5.1.4) 
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Table  5.1.1 

shows  the  results  from  one  computer  run  with 

the  generator 

x  -  5^^x 

i  +  1  "  5  xi 

(mod 

;  2&)  . 

Number  of 

Estimate  for  tt 

Tosses 

£  =  .011 

£  =  .005 

II 

• 

o 

VO 

a  =  .02 

a  =  .01 

a  =  .10 

10000 

3.1321 

3.1939 

3.1518 

20000 

3.1464 

3.1615 

3.1422 

30000 

3.1384 

3.1695 

3.1522 

40000 

3.1208 

3.1638 

3.1564 

50000 

3.1219 

3. 1693 

3.1565 

60000 

3.1260 

3.1585 

3.1571 

70000 

3.1191 

3.1598 

3.1517 

80000 

3.1161 

3.1535 

3.1535 

90000 

3.1185 

3.1502 

3.1515 

100000 

3.1215 

3.1592 

3.1521 

Table  5.1 

. 1  -  Estimates 

for 

it  Using  the  Buffon 

Needle 

Technique 

From  Table  5. 

1.1,  it  can  be 

seen 

that  only 

a  very  rough 

approximation 

to  tt  may  be  found 

using  the 

Buffon  approach 

If  it  were  possible  to  generate  infinitely  many  random 
numbers,  then  we  could,  theoretically,  calculate  tt  to  as 
many  significant  digits  as  we  desire.  However,  as  pointed 
out  in  Hull  and  Dobell  (1962),  the  solution  is  approached 


i  tf  fiuT  i9vjuqmoo  r  fto  rno' 
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with  an  error  of  order  n2  ,  where  n  is  the  number  of 

trials.  This  means  that  an  additional  significant  digit 
2 

requires  n  as  many  computations  as  the  previous  digit. 
It  can  be  seen  from  Table  5.1.1  that  only  two  significant 
digits  result  from  100,000  tosses  of  the  needle. 

5 . 2  Chuck-a-Luck 

One  example  of  a  game  of  chance  is  the  familiar 
chuck-a-luck .  The  game  is  played  with  three  ordinary 
dice  which  are  tossed.  A  player  is  paid  according  to  the 
number  of  dice  which  show  his  bet.  Table  5.2.1  shows  the 
probabilities  and  the  return  to  a  player  who  has  bet  one 
dollar . 


Result  of  Throw  Probability  Payoff 


No  6  ’  s 

.579 

-$1 

One  6 

.3^7 

$1 

Two  6 1 s 

.069 

$2 

Three  6fs 

.005 

$3 

1.000  -$  .  0787 

Table  5*2.1  -  Payoff  to  Player  Betting  on  6 

From  Table  5.2.1,  it  can  be  seen  that  a  player  will  expect 
to  lose  7.87  cents  whenever  the  dice  are  tossed.  From 
this,  we  can  show  that  a  player  who  starts  with  x  dollars 
can  expect  to  have  lost  all  his  money  after 
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N  = 


x 

.  07B7 


(5.2.1) 


throws  of  the  three  dice.  The  expected  duration  of  the 
game  is  shown  in  Table  5.2.2. 


Initial 

Expected 

Mimimum 

Maximum 

Amount 

Duration 

r 

Observed 

Observed 

$5 

63.5 

5 

385 

$10 

127.05 

12 

465 

$20 

254.1 

42 

2057 

Table  5.2.2 

-  Expected 

Duration  of  Chuck- 

a-Luck 

In  order  to  simulate  chuck-a-luck ,  four  random 
numbers  in  the  interval  (0,1)  were  generated.  The 
interval  was  divided  into  six  subintervals  to  determine 
the  player’s  bet  and  the  result  of  the  throw  of  the 
three  dice.  The  player's  payoff  was  determined  according 
to  the  number  of  dice  showing  his  bet.  The  results  for  a 
simulation  of  fifty  games  are  shown  in  the  Appendix. 

5 • 3  Integration  Under  a  Curve 

The  third  example  of  Monte  Carlo  techniques  is 

Id 

integration  under  a  curve,  i.e.,  /  f(x)dx  .  The  method 

a 

may  be  developed  in  the  following  manner. 


>- 
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The  maximum  value  of  the  function  f(x)  within  the  interval 
is  computed  and  given  a  value  Max  .  The  value  of  b-a 
is  computed  (where  a  and  b  are  the  lower  and  upper 
limits  of  integration  respectively)  and  given  a  value  of  T  . 
Two  random  coordinates  (x,y)  are  generated  within  the 
interval  (0,1)  and  transformed  as  follows: 


x  <-  x  Max 

y  y  T  . 

If  y  <_  f(x)  ,  then  the  total  number  of  points  on  or  below 
the  curve  is  incremented  by  one.  Otherwise,  the  total 
number  of  points  above  the  curve  is  incremented  by  one. 

The  ratio  of  the  value 


_ Total  points  on  or  below  the  curve _ 

Total  points  below,  on,  or  above  the  curve 


is  computed.  The  area  may  be  approximated  by 


Area  -  R  T  Max  . 
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The  number  of  points  needed  to  calculate  the  next 

2 

significant  digit  is  n  as  many  as  for  the  previous 
digit  -  a  situation  analogous  to  that  of  Buffon's 
needle . 

5 • 4  Hypothesis  Testing 

Specifically,  we  are  concerned  with  the  probability 
of  analyzing  a  disease  D  correctly  when  given  a  set  of 
symptoms  S  =  {  s^  ,  s^  ,  .  „  .  ,  s^.}  .  By  Bayesf  Theorem  (see 
Freund  (1962)),  we  know  that 


P  ( D  |  S  )  =  P  (  S  |  D )  . 


(5.4.1) 


Two  different  cases  could  be  examined  -  either  the 
symptoms  are  independent  of  each  other,  or  that  some  of  the 
symptoms  are  pairwise  dependent.  For  simplicity,  we  will 
consider  three  symptoms  s^,  s^,  and  associated 

estimates  of  the  conditional  probabilities  p^,  p2,  p^  of 
the  disease  D  given  the  symptoms  s^,  s^,  s^,  e.g„,  P^ 

is  an  estimate  of  p^  =  P(s1 |D)  and  P^  is  an  estimate 
of  p^^  =  P( s2ns^ | D)  . 

In  the  simulation,  let  us  consider  the  second  case, 
i o e „ ,  two  symptoms  2  and  3  are  pairwise  dependent.  Then, 


if  we  assume 


I 

>■ 


74 


where  is  the  frequency  of  the  symptom  s^  and  N 

the  number  of  observations  and 


f 

2n3 

23  N 


With  different  values  of  for  different  samples,  we 

should  be  able  to  obtain  the  distributions  of  P-^  and 
P^  •  We  shall  assume  that  p^  is  binomial  with  mean 
P^  and  variance 


Pi(1-V 

N 


and  p^^  is  binomial  with  mean  P^^  and  variance 


(1-P23) 


N 


is 


r  ■ 
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In  addition  to  obtaining  the  distributions,  we  would  like 
to  examine  the  results  of  two  different  probability 
estimates 


P  *^P 12^13^23 


(5.4.1) 


and 


P  3  ( P  2_P  2  3  2^  1 3  3^  1 2  ^  (5*4.2) 

of  detecting  the  disease  D  given  the  symptoms  s^3  s^, 

The  simulation  was  run  on  the  7040  computer  for  six 
different  conditions.  The  values  p^  and  p^  were  set 
constant  at  .5, however,  P3I2  was  varied  as  *5,  .7,  and  .9. 
The  sample  frequencies  of  PT  and  P"  were  calculated 
using  five  hundred  simulations  for  sample  sizes  of  twenty 
or  fifty  observations.  The  simulation  consisted  of 
generating  pseudo-random  numbers  with  a  uniform  distribution 
determining  whether  a  symptom  s-^,  s^,  or  s^  was  present 
and  tallying  the  frequencies.  The  results  are  shown  in 
Table  5.4.1. 


•.  ;  tj:  s 


tf98  e'isw 
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Sample  Size 

p  , 

3  |  2 

Pf 

pM 

20 

.5 

.1298 

.1295 

20 

.7 

.  1462 

.1409 

20 

.9 

.1664 

.1567 

50 

.5 

.  1268 

.1266 

50 

.7 

.1502 

.1443 

50 

.9 

.  1662 

.1569 

Table  5.4.1  - 

Results  from 

Simulation  of 

Disease 

Detection  for 

Symptoms  s^. 

s  2  }  s  ^  « 

The  simulation  may  be  carried  further  to  examine  binomial 
product  distributions  of  Pf  and  P"  for  the  entire 
set  of  500  observations. 

It  might  be  of  interest  to  examine  the  bias  in  the 
two  estimates  (5.4.1)  and  (5.4.2)  of  P-^g  *  ^  we 

assume  that  1  and  2  are  independent,  i.e.,  P-^  =  P-^2  5 

but  that,  in  general,  2  and  3,  are  not,  i.e., 


P  =  P  P 

r2  3 


2  3  |  2  • 
=  pnp„p. 


Then,  p’  and 


are  both  estimates  of 


123  1  2  3  |  2  * 

Substituting  population  values  in  (5.4.1),  we  get 


77 


/pi2p23pi3  y/(p1p2)  (p2P3  |2}  (^pip3) 


=  /(p1p2p3|2)2(p3/p3|2) 

P12 3  ’/p3/p3|2 


The  expression  under  the  radical  is  the  bias  which  is  1 
if  and  only  if  2  and  3  are  independent,  i.e.,  if  and  only 


Substituting  population  values  in  (5.4.2),  we  get 


3  ('P1P23+P2P13+P3P12  ^  “  3  ^  PlP2P3  |  2+P2PlP3+P3PlP2  ^ 


=  p-j_P2  ( P3  |  2+2P3 ) 


=  p 


123  v  3p 


P3|2+2P3 


3  |  2 


) 
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Again,  the  estimate  is  unbiased  if  and  only  if  ?3  =  P3|2 
Both  bias  terms  can  be  related  to  the  parameters 
of  the  simulation  by  substituting 


P 


3 


)  . 


P2  P3 I  2  +  (1_P2)  (1_P3|2 
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APPENDIX 


Two  multiplicative-congruent ial  pseudo-random  number 

generators  are  examined  in  the  Appendix.  The  first  uses  a 

modulus  of  2  ,  the  second  a  modulus  of  2-31  (the 

largest  prime  less  than  2  ) .  The  MAP  source  statements 

for  both  generators  are  shown.  Either  generator  may  be 

called  from  a  FORTRAN  mainline  program  in  the  same  manner 

as  a  FORTRAN  function  subprogram,  i.e.,  Y  =  RANDOM  (X)  , 

X  is  the  initial  value  to  be  used.  The  pseudo-random 

number  is  returned  as  a  real  constant  in  the  interval  (0,1). 

If  integer  pseudo-random  numbers  are  required  between  1  and 
35  35 

M  (,M=2  or  2  -31),  the  revisions  necessary  are  shown. 

The  frequencies  and  serial  frequencies  for  blocks 
of  10000  random  numbers  are  shown  in  the  next  section  of 
the  Appendix.  It  was  seen  that  both  methods  gave  short 
cycles  in  the  least  significant  digits,  so  only  nine  or 
ten  octal  digits  were  examined.  In  all  cases,  the  bit 
for  the  sign  was  ignored  for  the  reasons  stated  in  section 
2.5.  Both  generators,  with  the  appropriate  multipliers, 
gave  acceptable  chi-square  values. 

2 

The  results  from  Gorenstein's  test  and  the  d 
test  are  tabulated  in  the  next  section.  In  either  test, 
the  mult iplicat ive-congruent ial  method  gave  acceptable 
chi-square  values. 

Plots  of  the  Box  and  Muller  normal  transformation 


>■ 

' 
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and  Juncosa's  transformation  for  sixteen  and  sixty-four 
elements  are  given. 

The  last  section  of  the  Appendix  displays  the  results 
of  the  applications  discussed  in  Chapter  V. 


' 


.JlBMAP  NULCGN 
* 
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*  NL'LTIPIICATI  VE-CCNCRUENT  I  AL  PSEUDO  RANDOM  NUMBER  GENERATOR 

* 


*  ...  CALLING  SEQUENCE  FROM  FORTRAN  MAINLINE  IS  ... 

*  Y=R ANCON ( X ) 

*  WHERE  X  IS  THE  STARTING  VALUE 

*  CHOSEN  BY  THE  PROGRAMMER. 

*  IF  X  IS  EVEN,  IT  IS  MADE  ODD. 


* 

* 

* 


* 

ENTRY 

RANCCN  SAVE 
CLA 
T  PL 

$ 


THE  PSEUDO-RANDOM  NUMBER  RETURNED  HAS  A  VALUE  IN 

THE  OPEN  INTERVAL  (0,1). 


RANCCN 

1,2  ,  A 

SWITCH 
NOTH  ST 


DEFINE  THE  ENTRY  POINT 
SAVE  THE  RE  I  URN  LINKAGE 
TEST  FOR  FIRST  TIME 
THROUGH 


FIRST  TIME  THROUGH  GENERATOR  -  INITIALIZE  X 

MSP  SWITCH  RESET  SWITCH 

CLA*  3, A  GET  THE  INITIAL  VALUE  OF  X 

CRA  =CCCC0CCCC0C01  MAKE  X  ODD 

STO  X  STCRE  INITIAL  X 


* 

NCTFST 


* 

❖ 


END  OF  INITIALIZATION  -  MAIN  ROUTINE  FOLLOWS 


LCQ 

MPY 

STQ 

CLA 

ARS 

CRA 

FAD 

RETURN 


=1220703125 

X 

X 

8 

FLCT 
FLC  T 
RANCON 


PUT  X  INTO  MQ 
MULTIPLY  BY  5**13 

REPLACE  OLD  X  WITH  NEW  X  (MOD  2**35) 
PUT  X  INTO  AC 

*SE  T  UP  X  FOR  FLOATING  POINT 
*  INTRODUCE  EXPONENT 
^NORMALIZE  MANTISSA 
RETURN  TO  MAINLINE  WITH  X  IN  AC 


CONSTANTS 


❖ 


X  BSS 
FLOT  CCT 
SW  I  TCH  CCT 


l  STORAGE  FOR  X 

‘2CCCCCC0CCCC  FLOATING  POINT  EXPONENT  /  NORMAL  I  Z  E  R 
ACCCCCCCCCCO  INITIALIZE  AS  A  NEGATIVE  NUMBER 


*  INSTRUCTION  DESCRIPTIONS  BEGINNING  WITH  *  MAY  BE  DELETED 

*  IF  INTEGER  VALUES  BETWEEN  l  AND  2**35  ARE  REQUIRED 

*  INSTEAD  OF  FLOATING  POINT  NUMBERS 


END 


END  OF  SUBROUTINE 


-  rr  1  3 

...  =  5  x. 
l+l  i 


(mod  235) 


MAP  Listing  for  x 


i  I  BMAP 
* 

* 

* 

❖ 

* 

* 

* 

* 

* 

* 

RANLEH 


❖ 

* 

❖ 


LEMPER 
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* 

* 


PULTIPLICAT I VE-CONGRUENT  IAL  PSEUDO-RANDOM  NUMBER  GENERATOR 
USING  A  MODULES  CF  THE  LARGEST  PRIME  PCSSIBLE 
IN  THE  7C4G 

...  CALLING  SEQUENCE  FRCP  FORTRAN  MAINLINE  IS  ... 

Y  =  R  ANL EH { X  ) 

VALUE  RETURNED  IS  SAME  AS  FOR  MODULUS  2**35  PROGRAM 


ENTRY 

SAVE 

CLA 

TPL 


RANLEH 
1,2, A 
SW  I  TCI! 
NCTFST 


INITIALIZATION  FCR  X 


PSP 
C  L  A: 
S  TO 


SWITCH 
3,  A 
X 


DEFINE  THE  ENTRY  POINT 
SAVE  THE  RETURN  LINKAGE 
IS  THIS  THE  FIRST  TIME 
NC,  TRANSFER  TO  NOTES T 


RESET  SWITCH 
GET  THE  INITIAL  X 
AND  STORE  IT 


❖ 

❖ 

* 

NCTFST 


ROUTINE  TC  CALCULATE  NEW  X 

LDQ  X  PUT  X  INTO  MQ 

PPY  =1220703125  MULTIPLY  BY  5**13 

CVP  =03 /7 7 77777741  DIVIDE  TO  GET  MOO  2**35-31 

STO  X  STORE  NEW  X 

ARS  8  *  S  F  I F  T  X  TO  FLOAT  IT 

ORA  FLCT  *  INSERT  EXPONENT 

FAD  FLCT  ^NORMALIZE  IT 

FAD  FLCT  ^NORMALIZE  THE  PSEUDO-RANDOM  NUMBER 

RETURN  RANLEH  RETURN  TO  MAINLINE 

CONSTANTS 


X 

FLOT 

SWITCH 


ess 

CCT 

ccr 


l 

2CCCC0CCCCC0 

4CCCCCCCCCCC 


STORAGE  FOR  X 

FLOATER  AND  NORMAL  I ZER 

SWITCH  FCR  FIRST  TIME  THROUGH 


❖ 

* 


INSTRUCTIONS  WITH  *  IN  CC30  MAY  BE  OMITTED  IF  AN  INTEGER 
BETWEEN  1  AND  2**35-31  IS  DESIRED 


END 


END  OF  SUBROUTINE 


=  5‘^x.  (mod  235_3i) 
l  +  l  l 


MAP  Listing  for  x 


>- 
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serial  FREQUENCY  FOR  octal  POSITION  I 


0 

1 

2 

3 

4 

5 

6 

7 

c 

148 

149 

164 

153 

147 

1  79 

153 

165 

1 

152 

l  3d 

156 

155 

16  1 

1  76 

159 

141 

2 

149 

168 

172 

147 

152 

16  7 

159 

147 

3 

160 

170 

150 

170 

129 

130 

189 

1  72 

4 

164 

14  1 

150 

157 

147 

144 

145 

143 

5 

149 

16  1 

147 

170 

156 

150 

164 

16  3 

6 

167 

158 

160 

1  70 

161 

168 

153 

142 

7 

14  9 

153 

162 

168 

1  38 

146 

157 

150 

1258 

1238 

126  l 

1290 

1191 

1260 

1279 

1223 

CF  I-SGL ARE  FOR 

SERIAL 

FREQUENCY  IS 

58.52 

CFI- 

SQUARE  FOR 

FREQUENCY  IS 

5.66 

SERIAL  FREQUENCY  FOR  OCTAL  POSITION  2 


0 

1 

2 

3 

4 

5 

6 

7 

C 

163 

167 

157 

147 

164 

166 

1  72 

150 

1 

170 

14  7 

159 

150 

14  7 

139 

150 

153 

2 

154 

143 

14? 

158 

l  74 

162 

155 

170 

3 

161 

146 

154 

163 

134 

164 

160 

134 

4 

150 

1  59 

17  1 

145 

190 

138 

151 

174 

5 

16  6 

144 

166 

152 

150 

142 

147 

169 

6 

167 

164 

15  1 

149 

144 

166 

153 

158 

7 

155 

145 

158 

152 

175 

159 

164 

15  1 

1286 

12  15 

1253 

1216 

1278 

1236 

12  52 

1259 

cfi-scuare  FOR 

SERIAL 

FREQUENCY  IS 

49.7  7 

CFI- 

SGUARE  FOR 

FREQUENCY  IS 

3 . 85 

SERIAL  FREQUENCY 

FOR  OCTAL 

POSITION  3 

0 

l 

2 

3 

4 

5 

6 

7 

C 

147 

153 

151 

14  8 

152 

146 

15  5 

l  76 

l 

154 

176 

154 

160 

169 

183 

170 

151 

2 

145 

180 

154 

162 

157 

149 

141 

162 

3 

141 

166 

164 

125 

140 

128 

16  3 

165 

4 

160 

144 

16  7 

154 

143 

144 

153 

158 

5 

104 

155 

156 

149 

150 

149 

151 

149 

6 

169 

17  3 

134 

132 

l  5  4 

169 

169 

157 

7 

128 

170 

170 

162 

158 

1  75 

155 

172 

1228 

13  17 

1250 

1  192 

1223 

1243 

1257 

1290 

CF  I-SGL ARE  FOR  SERIAL  FREQUENCY  IS  71.05 
CHI-SCCARE  FOR  FREQUENCY  IS  8.61 


Serial  Frequency  Results  for  xj_  +  i  E  5'  (mod  2  ) 


O  VJ1  ^  IjJ 
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SERIAL  FREQUENCY  FOR  OCTAL  POSITION  A 


C 

1 

2 

3 

4 

5 

6 

7 

c 

139 

150 

151 

145 

181 

139 

157 

168 

1 

165 

162 

164 

144 

155 

1  o9 

1  60 

156 

2 

141 

158 

154 

144 

163 

1  77 

158 

139 

3 

164 

148 

147 

152 

14  5 

166 

128 

166 

A 

140 

153 

148 

137 

149 

171 

152 

185 

5 

168 

169 

150 

169 

156 

148 

160 

170 

f : 

148 

168 

160 

150 

137 

169 

166 

136 

7 

165 

14  7 

160 

175 

149 

171 

153 

146 

l  230 

12  7  5 

1234 

12  16 

1235 

1310 

1234 

1266 

CH1-SGLARE  FOR 

SERIAL 

FREQUENCY  IS 

6  5.75 

CF  I- 

SQUARE  FOR 

FREQUENCY  IS 

5.42 

SERIAL  FREQUENCY  FOR  OCTAL  POSITION  5 


0 

1 

2 

3 

4 

5 

6 

7 

C 

163 

15  1 

172 

145 

149 

151 

169 

143 

1 

146 

149 

148 

154 

150 

164 

180 

156 

2 

161 

153 

17  7 

14  4 

15  2 

153 

146 

172 

3 

159 

160 

152 

157 

161 

155 

160 

159 

4 

147 

155 

147 

168 

143 

151 

163 

146 

5 

160 

1  74 

171 

164 

153 

164 

146 

137 

6 

153 

159 

147 

169 

14  7 

162 

161 

169 

7 

154 

146 

144 

162 

165 

169 

142 

151 

1243 

12  4  / 

1258 

1263 

1220 

1269 

1267 

1233 

CFI-SCUARE  FOR 

SERIAL 

FREQUENCY  IS 

38.35 

CHI- 

SGUARE  FOR 

FREQUE 

NCY  I  S 

1 . 70 

SERIAL  FRECLENCY  FOR  OCTAL  POSITION  6 


0 

1 

2 

3 

4 

5 

6 

7 

146 

153 

135 

159 

169 

164 

157 

143 

124 

136 

165 

176 

15  5 

15  5 

155 

147 

139 

142 

160 

163 

1  54 

156 

147 

161 

156 

148 

169 

146 

170 

149 

156 

164 

171 

165 

148 

168 

180 

165 

185 

158 

173 

17  1 

142 

136 

17  7 

152 

168 

136 

177 

14  l 

162 

150 

180 

160 

151 

150 

140 

157 

14  1 

160 

15  5 

154 

152 

156 

1226 

12  13 

1222 

1258 

1340 

1255 

1271 

12  15 

CHI-SQUARE  FOR 

SERIAL 

FREQUENCY  IS 

65.83 

CHI-SQUARE  FOR  FREQUENCY  IS  10.07 
Serial  Frequency  Results  for  xi+q  E  5  '“x^Cmod  2  ) 


c 

1 
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7 

C 
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6 
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SERIAL  FREQUENCY 

FCR  OCTAL 

POSIT IGN  7 

c 

1 

2 

3 

4 

5 

6 

7 

140 

168 

151 

151 

155 

148 

164 

148 

170 

157 

158 

163 

173 

1  74 

134 

161 

156 

163 

167 

137 

1  86 

162 

170 

155 

139 

168 

166 

158 

145 

146 

154 

133 

164 

151 

171 

154 

152 

149 

155 

177 

168 

157 

167 

154 

140 

168 

143 

158 

152 

168 

162 

140 

163 

152 

166 

1  32 

136 

158 

154 

152 

159 

156 

149 

153 

1225 

1250 

1296 

1209 

1273 

1255 

1235 

1217 

CF I-SGU ARE  FOR  SERIAL  FREQUENCY  IS  53.24 
CHI-SCUARE  FCR  FREQUENCY  IS  6.31 


SERIAL  FREQUENCY 

FCR  OCTAL 

POSITION  8 

0 

1 

2 

3 

4 

5 

6 

7 

152 

155 

168 

142 

148 

163 

148 

169 

149 

152 

168 

157 

175 

141 

156 

178 

153 

167 

153 

146 

167 

147 

149 

161 

147 

168 

150 

148 

165 

146 

175 

145 

145 

170 

150 

168 

146 

144 

163 

158 

161 

156 

142 

164 

147 

150 

167 

147 

170 

148 

148 

168 

149 

173 

144 

152 

163 

160 

164 

151 

147 

l  70 

150 

152 

1245 

1276 

1243 

1244 

1244 

1234 

1252 

1262 

CFI-SQUARE  FCR  SERIAL  FREQUENCY  IS  41.11 
CF  I -  SQUARE  FCR  FREQUENCY  IS  0.98 


SERIAL  FREQUENCY 

FOR  OCTAL 

PUS  If  ION  9 

0 

I 

2 

3 

4 

5 

6 

7 

159 

156 

177 

158 

137 

175 

138 

157 

157 

136 

176 

136 

155 

157 

156 

174 

139 

155 

158 

156 

175 

158 

1  36 

178 

15-6 

176 

157 

136 

176 

137 

155 

156 

136 

176 

136 

157 

155 

157 

174 

154 

158 

155 

157 

177 

157 

138 

175 

1  37 

176 

156 

137 

173 

134 

157 

155 

156 

176 

137 

15  7 

156 

156 

175 

155 

137 

1257 

1247 

1255 

1249 

1245 

12  54 

1244 

1249 

CFI-SQUmRE  FCR  SERIAL  FREQUENCY  IS  77.93 
CHI-SQUARE  FCR  FREQUENCY  IS  12.96 

1 R  35 

Serial  Frequency  Results  for  xj_  +  q  E  5  Jx_^  (mod  2  ) 


- 
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TOTAL  SERIAL  FREQUENCY 

0 

l 

2 

3 

A 

5 

6 

7 

L  35  l 

1  AC  2 

1  A  2  6 

1  3  A  8 

1A02 

1  A3  1 

1413 

1419 

1387 

1353 

1 A  A  8 

1395 

1 A  A  0 

1  A58 

1 A  2  0 

1417 

1337 

l  A  2  9 

1 A  3  7 

1357 

1 A  8  0 

l  A3  l 

1361 

l  A  A  5 

1  A  0  3 

1 A  5  0 

1  AC  9 

1355 

1  3 1>  5 

1321 

1  A  AO 

1394 

1377 

1 A  1  A 

1388 

lAoa 

1  AO  5 

1  363 

1 A  A  1 

14  5  3 

1  A  8  7 

1 A  6  2 

1398 

1  A3  5 

1386 

1  361 

1  A2  1 

1366 

1 A  79 

1 A  3  5 

1361 

1 A  0  1 

1369 

1  A  76 

1418 

1352 

J  371 

13  73 

1A  10 

l  A3  8 

1 A02 

1 A  7  5 

1377 

1368 

1  1198 

113  18 

1  12  77 

1  1  L  3  7 

1 1 2  A  9 

11316 

11291 

11214 

CF  I-SGL ARE  FOR 

SERIAL 

FREQUENCY  IS 

70. 16 

CHI-SQUARE  FCR  FREQUENCY  IS  2.50 


x 


i  +  1 


Serial  Frequency  Results  for 


x . 
1 


(mod  235) 


c 

1 

2 

3 

4 

5 

6 

7 

C 
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2 

3 

4 

5 

6 

7 

C 

l 

2 

3 

4 

5 

6 
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SERIAL  FREQUENCY  FCR  OCTAL  POSITION  L 


0 

l 

2 

3 

4 

5 

6 

7 

139 

146 

147 

170 

167 

154 

1  36 

153 

172 

152 

144 

179 

159 

165 

155 

133 

139 

172 

170 

151 

160 

139 

16  7 

143 

174 

177 

165 

165 

164 

190 

160 

145 

148 

146 

162 

164 

157 

161 

165 

163 

16  1 

147 

148 

192 

157 

141 

150 

166 

142 

146 

154 

156 

162 

172 

144 

149 

137 

173 

151 

163 

140 

140 

148 

143 

1212 

1259 

124  l 

l  340 

1266 

1262 

1225 

1195 

CFI-SCLARE  FCR  SERIAL  FREQUENCY  IS  69.81 
CF I -SQUARE  FCR  FREQUENCY  IS  11.00 


SERIAL  FREQUENCY 

FCR  OCTAL 

POSITION  2 

0 

1 

2 

3 

4 

5 

6 

7 

161 

149 

16  6 

140 

149 

149 

166 

158 

151 

138 

165 

178 

163 

147 

133 

151 

152 

155 

160 

171 

152 

150 

160 

173 

181 

159 

163 

L  59 

171 

148 

165 

143 

161 

162 

133 

150 

l  50 

144 

164 

169 

14  1 

150 

189 

153 

1  36 

154 

164 

125 

128 

159 

134 

162 

166 

167 

169 

180 

163 

154 

163 

1  76 

146 

153 

144 

165 

1238 

1226 

1273 

1289 

1233 

1212 

1265 

12  64 

CFI-SCLARE  FCR  SERIAL  FREQUENCY  IS  70.49 
CFI-SQUARE  FCR  FREQUENCY  IS  3.94 


SERIAL  FREQUENCY  FOR  OCTAL  POSITION  3 


0 

l 

2 

3 

4 

5 

6 

7 

159 

18  1 

141 

154 

156 

157 

164 

173 

177 

2C  1 

167 

180 

168 

165 

150 

152 

144 

169 

172 

154 

156 

157 

152 

152 

153 

173 

156 

150 

156 

177 

140 

146 

156 

165 

158 

147 

149 

151 

14  7 

152 

153 

153 

165 

153 

165 

173 

149 

143 

163 

155 

148 

157 

127 

144 

135 

148 

180 

163 

149 

156 

148 

130 

140 

126 

1285 

1360 

1256 

125  1 

1225 

12  54 

1177 

1192 

CFI-SQUARE  FCR 

SERIAL 

FREQUENCY  IS 

74.23 

CHI-SQUARE  FCR  FREQUENCY  IS  18.17  * 

13  35 

Serial  Frequency  Results  for  xj_  +  q  E  5  (mod  2  -31) 


. 

c 

l 

2 

3 

4 

5 

6 

7 

C 

1 

2 

3 

4 

5 

6 

7 

C 

1 

2 

3 
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5 

6 
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SERIAL  FRFGUENCY  FOR  OCTAL  POSITION  4 


0 

l 

2 

3 

4 

5 

6 

7 

147 

180 

13  7 

148 

133 

153 

136 

149 

151 

150 

144 

142 

160 

164 

172 

174 

127 

155 

144 

152 

155 

158 

143 

173 

141 

155 

153 

141 

191 

148 

176 

159 

157 

136 

143 

199 

167 

138 

163 

167 

177 

16  0 

149 

153 

155 

161 

133 

158 

129 

149 

156 

165 

150 

162 

140 

1  78 

159 

172 

176 

164 

159 

162 

16b 

176 

1  188 

125  7 

1207 

1  264 

12  7  5 

124  6 

1229 

1334 

CP  I-SQLARE  FOR  SERIAL  FREQUENCY  IS  86.35  * 
CHI-SQUARE  FOP  FREQUENCY  IS  11.26 


SERIAL  FREQUENCY  FOR  OCTAL  POSITION  5 


0 

1 

2 

3 

4 

5 

6 

7 

156 

137 

166 

16  3 

165 

149 

138 

157 

154 

164 

147 

178 

1  54 

151 

160 

162 

151 

170 

138 

152 

158 

150 

165 

156 

158 

17  7 

168 

162 

153 

161 

162 

157 

161 

152 

164 

173 

156 

133 

123 

178 

166 

151 

162 

172 

152 

164 

144 

136 

143 

160 

143 

140 

151 

170 

150 

143 

142 

159 

152 

163 

151 

169 

153 

175 

1231 

1 2  7  C 

1240 

1303 

124  5 

1247 

1200 

1264 

CHI-SQUARE  FOR 

SERIAL 

FREQUENCY  IS 

52.  16 

CH  I-SQ 

U ARE  FCR 

FREQUENCY  IS 

5.  12 

SERIAL  FREQUENCY 

FOR  OCTAL 

POSITION  6 

0 

l 

2 

3 

4 

5 

6 

7 

168 

192 

160 

156 

159 

156 

181 

136 

150 

155 

171 

158 

154 

144 

163 

139 

164 

137 

166 

159 

14  7 

14  1 

182 

170 

160 

156 

157 

165 

1  72 

128 

158 

152 

157 

157 

157 

160 

163 

154 

141 

157 

182 

140 

147 

148 

14  1 

147 

151 

143 

16  8 

14  8 

161 

141 

176 

169 

176 

156 

159 

149 

147 

161 

128 

166 

143 

145 

1308 

1234 

1266 

1248 

1246 

1  205 

1295 

1198 

CHI-SQUARE  FOR  SERIAL  FREQUENCY  IS  67.87 
CHI-SQUARE  FOR  FREQUENCY  IS  3.52 

1 R  35 

Serial  Frequency  Results  for  x. ,,  =  5  x.  (mod  2  -31) 

i+j-  l 
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SERIAL  FREQUENCY 

FOR  OCTAL 

POSITION  7 

0 

l 

2 

3 

4 

5 

6 

7 

146 

169 

154 

155 

147 

157 

160 

163 

156 

143 

169 

161 

177 

142 

152 

144 

155 

157 

129 

152 

168 

154 

146 

165 

161 

165 

174 

162 

156 

143 

154 

155 

167 

160 

165 

154 

157 

168 

169 

145 

146 

156 

133 

169 

147 

145 

157 

178 

167 

160 

161 

136 

1  76 

160 

158 

134 

153 

134 

141 

161 

157 

162 

156 

157 

1251 

12  44 

1226 

1270 

1285 

1231 

1252 

1241 

CHI-SQUARE  FOR  SERIAL  FREQUENCY  IS  53. CO 
CFI-SCUARE  FOR  FREQUENCY  IS  2.15 


SERIAL  FREQUEtNCY 

FOR  OCTAL 

POSITION  8 

0 

1 

2 

3 

4 

5 

6 

7 

178 

165 

150 

143 

163 

155 

150 

155 

159 

176 

173 

164 

163 

160 

154 

181 

15  7 

185 

16  4 

155 

l  70 

155 

112 

158 

165 

152 

148 

152 

148 

136 

158 

142 

155 

18  1 

163 

153 

171 

158 

155 

154 

143 

148 

140 

164 

160 

144 

156 

141 

140 

160 

167 

128 

156 

144 

155 

158 

162 

16  3 

151 

142 

159 

144 

168 

l  71 

1259 

1 3  3  C 

1256 

1201 

129  0 

1196 

1208 

12  60 

CF I~ SQUARE  FOR  SERIAL  FREQUENCY  IS  65.31 
CHI-SQUARE  FOR  FREQUENCY  IS  L2.23 


SERIAL  FREQUENCY 

FOR  OCTAL 

POSITION  9 

0 

l 

2 

3 

4 

5 

6 

7 

174 

164 

167 

154 

159 

145 

149 

157 

133 

168 

1  5  1 

147 

176 

164 

14  1 

167 

158 

1  7  7 

159 

155 

157 

138 

141 

185 

156 

150 

175 

154 

135 

167 

139 

150 

150 

154 

146 

136 

133 

152 

165 

186 

161 

142 

14  7 

169 

15  1 

166 

162 

149 

163 

133 

162 

143 

156 

149 

163 

143 

174 

159 

163 

168 

155 

1  6  6 

152 

170 

1269 

124  7 

1270 

1226 

1222 

1247 

1212 

1307 

CHI-SQUARE  FOR  SERIAL  FREQUENCY  IS  64. 3L 
CHI-SQUARE  FCR  FREQUENCY  IS  5.47 


1  o  u 

Serial  Frequency  Results  for  x^  +  ]_  5  5  Jx^  (mod  2  -31) 
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SERIAL  FREQUENCY  FOR  OCTAL  POSITION  10 


0 

l 

2 

3 

4 

5 

6 

7 

173 

150 

165 

177 

17  1 

14  6 

151 

171 

147 

143 

177 

151 

162 

145 

156 

172 

171 

1  7  4 

159 

155 

169 

136 

156 

142 

179 

148 

158 

168 

145 

158 

169 

147 

136 

17  1 

167 

152 

161 

146 

165 

161 

154 

127 

157 

154 

141 

144 

146 

151 

177 

160 

L  3  3 

183 

1  50 

150 

1  32 

143 

167 

180 

146 

132 

160 

149 

153 

161 

1304 

12  5  3 

1262 

12  72 

12  59 

l  174 

1228 

124  8 

CF  I-SQU ARE  FOR 

SERIAL 

FREQUENCY  IS 

7  2.50 

CF  I- 

SQUARE  FOR 

FREQUENCY  IS 

7.92 

TOTAL  SERIAL  FREQUENCY 

0 

1 

2 

3 

4 

5 

6 

7 

1601 

1633 

1553 

1560 

15  74 

1521 

1531 

L  5  7  2 

1550 

1590 

1608 

1638 

16  3  6 

154  7 

1536 

L  5  7  5 

1518 

165  1 

156  1 

1556 

15  9  2 

14  7  8 

1524 

16  17 

1628 

16  12 

1617 

1578 

15  96 

1556 

1581 

1496 

1543 

1584 

1563 

1588 

1564 

1505 

1562 

1632 

15  84 

14  7  4 

153  7 

1627 

15  1  l 

1  539 

1512 

1490 

1520 

1530 

1519 

151  1 

1570 

158  7 

1522 

15  32 

1596 

1 6  C  6 

1539 

1606 

1503 

1541 

1  523 

1589 

12  545 

12680 

1249  7 

12  664 

12546 

12274 

12291 

12503 

CF  1-SQUARE  FOR 

SERIAL 

FREQUENCY  IS 

76.04 

CHI-SCUARE  FOR  FREQUENCY  IS  12.65 

1 R  RE 

Serial  Frequency  Results  for  xi+1  E  5  (mod  2-31) 
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AVERAGE 

2ND 

3RC 

TOTAL 

DIFF  IN 

LONGEST 

MOMENT 

MOMENT 

RUNS 

S.D. • S 

RUN 

C  .  5  C  3  4  16 

0.33672 

C.253C8 

6659 

0. 17394 

6 

0.5CC746 

0.33310 

C. 2 49  19 

6  716 

1  .  1  7806 

6 

0.492954 

0.32683 

C. 2 44 4 6 

6674 

0.18185 

6 

0. 50 123 3 

0.33441 

0.25115 

6  6  64 

0.05534 

6 

0 . 5C8 156 

0 . 3  4  C  0  7 

C  .255  17 

66  75 

0.20557 

6 

0.502428 

C  .  3  3  53  1 

0.25121 

6708 

C. 98830 

6 

0.496034 

G.  32891 

0.246C1 

6706 

C. 94086 

6 

0.497888 

0 . 3  3  2  8  C 

C. 25062 

6637 

0.69576 

6 

0.5C2285 

0.33641 

C  .25367 

6610 

1 .336  18 

6 

0.5C4452 

0.33802 

C.254C6 

6712 

1.083  18 

7 

0.495C96 

0.33263 

0.24982 

66  79 

0.30045 

6 

0 .5C1351 

0.33424 

C. 25049 

67  17 

1.20178 

6 

0.499638 

0.33215 

0.24857 

6724 

1.36781 

6 

0.499202 

0.33358 

C. 25061 

6  640 

0.62460 

6 

0.497473 

0.33025 

0.24677 

6  b  5  7 

0.221 38 

6 

0.5CC909 

0.33333 

0.24939 

6  7  11 

1.05946 

6 

0.493650 

C.3271C 

C  .24449 

6638 

0.6 7204 

7 

0.499264 

0.33237 

C  .24901 

6741 

1.  7  7104 

8 

0.499202 

0.33121 

0.247C6 

6642 

0.57717 

6 

0.497345 

0.33114 

0.24846 

6699 

0.77483 

7 

0 . 50 1354 

0.33494 

0.25171 

6  6  67 

0.01581 

7 

0.499121 

0.33328 

C.25C4C 

6  6  61 

0.12650 

7 

0.497  7  17 

0 .33058 

0.247  1  1 

6694 

0.65623 

7 

0.5CC846 

0.33472 

0.2515C 

6654 

0.29254 

6 

0.498689 

C  .  3  3  1 4  5 

0.24786 

6606 

1 .43106 

9 

0 .502060 

0.33583 

0.25227 

6695 

0.67995 

6 

0.501654 

0.33376 

0.24968 

6707 

0.96458 

7 

0.498795 

0.33277 

0.24996 

6642 

0.57717 

6 

0.5CC907 

0.33481 

C. 25141 

6  669 

0.06325 

7 

0.501666 

0.33444 

0.25088 

6701 

0. 82227 

6 

0 . 502848 

0.33622 

0.25284 

6593 

1.73941 

6 

0.497 779 

0 .33205 

C  .  2  4  9  3  0 

66  12 

1.28b 74 

6 

0.499485 

0 . 3  3  1 9  C 

C.248C5 

6  680 

0.324  16 

6 

C. 497324 

0 . 3  3  1 C  6 

C. 2 48  17 

6577 

2. 1 1892 

7 

0.501334 

0.33561 

0.25235 

6687 

0.49020 

6 

0. 504 8 82 

0.3374C 

0.25335 

6655 

0.26882 

6 

0.502317 

0.33608 

C  .25274 

6683 

0.39532 

6 

C. 5018  19 

0.33464 

0.25083 

6639 

0.64832 

8 

0.5C2347 

0.33696 

0.25384 

6753 

2.05567 

£ 

7 

0. 502  190 

0.33452 

0 .25002 

6640 

0.62460 

8 

0.508719 

0.34269 

0.25873 

6718 

1.22550 

- 

6 

0.501636 

0.33476 

C. 25101 

6629 

0.80552 

7 

0.496322 

0. 33C03 

0.24726 

6  7  50 

1.98451 

$ 

6 

0.5C2888 

0.33514 

0.25097 

6611 

1.31246 

5 

0.495526 

0.32878 

0.24589 

6  719 

1.24921 

7 

0 .503060 

0.33681 

C. 2 53 31 

6  6  86 

0.466^8 

6 

0.501794 

0.33577 

C. 25237 

6  691 

0.58508 

6 

0.498860 

0.33251 

0.24932 

6617 

1.17015 

7 

0.499773 

0. 33333 

0.25022 

6619 

1.12271 

7 
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Gorenstein's  Results  for  xi  +  ]_  5  5  (mod  2  ) 


- 

' 

98 


AVERAGE 

2N0 

3RC 

MCNEN  T 

MOMENT 

0.499439 

0.33339 

C  .25035 

0 . 4  99  7  30 

0.33365 

0.25096 

0. 492 30 8 

0 .32415 

C.24C67 

C. 496 851 

0.32955 

0.2 46 20 

0 .5C3C27 

0.33863 

C .25285 

0 .496214 

0.33013 

0.24730 

C . 5C4  1  2  2 

0.33811 

0.25497 

0.498  175 

0.33172 

0.24852 

0.49477G 

0 .328  10 

0.24670 

0 .499750 

0.33388 

C  .25082 

0 .501712 

0 .33490 

0.2512C 

0.498852 

C. 33198 

0.24872 

0.498639 

0.33257 

C  .24972 

0.494429 

0.32657 

0.24615 

0 . 5C3  30  G 

0.33779 

C. 25435 

0.502125 

0.33548 

0 . 2  5  2  C  3 

0.498457 

0.33252 

0.24983 

0.498404 

0. 3  3  20  5 

C  .249  12 

C  .499022 

0.33128 

0.24740 

0.502535 

0.33651 

0.25336 

0 . 503  102 

0.33627 

C  .25284 

0.501612 

0.3349  1 

C  .  2  5  1  1  9 

0.495857 

0 .32937 

0.24645 

0.494995 

0.32821 

0.24516 

0.49597C 

0. 32888 

C  .24566 

0.502487 

0.33492 

C  .25036 

C  .502  189 

C.  336C  / 

0.25309 

C. 5 0069 9 

0.33363 

C.  2 500  5 

0. 50C387 

0.33353 

0.25021 

0.496065 

0 .32933 

C. 24622 

0.499594 

0.33351 

0.25046 

0.488094 

0.33247 

C. 24962 

0 . 500  5  39 

0.33312 

0.24908 

0.502230 

0.33567 

0.25244 

0.492452 

0.32612 

0.24348 

0.498236 

0.33207 

0.24925 

0.503397 

0.33667 

C.253C0 

0.488036 

0 .33082 

0.24731 

C. 5006  1  7 

0.33514 

0.25208 

0 . 50 1 164 

0.33454 

C.2508  1 

0.501284 

0.33451 

0.25136 

0.502973 

0.33583 

C  .  25202 

0 .488702 

0.33133 

0.24762 

0. 50C668 

0.33336 

0.24977 

0.497858 

0.33086 

0.24751 

0.498660 

0.33275 

0.24982 

0.495982 

0 .32883 

0.24577 

0 . 504  8  76 

0 .33854 

0.25493 

0.496975 

0 .33031 

C  .24720 

Gorenstein’s  Results  for  x. 


TOTAL 

DIEE  IN 

LONGEST 

RUNS 

S.O. » S 

RUN 

6699 

0.77483 

6 

6697 

0.72739 

6 

6623 

1.02783 

6 

6  6  98 

0.75111 

6 

6619 

1.12271 

6 

6  7  15 

1.15434 

6 

6635 

0.74320 

6 

66  15 

1.21759 

6 

6  6  77 

0.25301 

6 

6684 

0. 4 1904 

8 

6623 

1 .02783 

6 

6  6  89 

0.53764 

7 

6636 

0.71948 

7 

6696 

0.70367 

6 

6  594 

1.71569 

6 

66  75 

0.20557 

6 

6699 

0.  / 7483 

6 

6687 

0.49020 

6 

6670 

0.08687 

7 

6694 

0.65623 

6 

6  6  74 

0.18185 

6 

6610 

1.33618 

7 

665  7 

0.22138 

7 

6762 

2.26914 

'r 

6 

6  768 

2.41146 

❖ 

6 

6653 

0.31625 

7 

6  70  7 

0.96458 

7 

6  748 

1.93707 

6 

6731 

1.53385 

7 

6675 

0.20557 

5 

6573 

2.21379 

'r 

5 

6655 

0.26882 

6 

6690 

0.56136 

7 

6  6  34 

0.76692 

6 

6720 

1.27293 

6 

66  79 

0.30045 

7 

6669 

0.06325 

6 

6623 

1.02733 

6 

6676 

0.22929 

6 

6642 

0.5/717 

7 

6654 

0.29254 

8 

6702 

0.84599 

8 

6683 

0.39532 

6 

6650 

0.38741 

7 

6715 

l . 1 5434 

6 

6  7  09 

1 .0 1202 

6 

6710 

1.03574 

6 

6681 

0.34788 

7 

6680 

0.324  16 

6 

-  6  x . 

(mod  235- 

-3D 

d2 

Expected 

Observed 

Accumulated 

Frequency 

Frequency 

Frequency 

C.  1C 

0.236332 

0. 231560 

U.  23  1560 

0 .20 

0.176973 

0.  1  75000 

0.607360 

0. 3C 

C  .  1396  95 

0.  139000 

0. 566360 

0.6C 

0,112/18 

0.  1  09160 

0.655520 

0 . 5C 

0 .090969 

0.091660 

0.76  71  60 

0.60 

0.072616 

0. 076880 

0.822060 

0.7C 

0.056769 

0. 056600 

0 . 878660 

0. 80 

0.062813 

0. 062080 

0.920520 

0 .90 

0.030631 

0. C32280 

0.952800 

1 .00 

0.019333 

0.019680 

0 . 9  7  26  80 

1  .  10 

C  .  0  1 0  7  7  7 

0.011920 

0.986600 

1  .20 

0. CCb365 

0.007660 

0.991080 

l  .  30 

0 . 0  C  3  7  6  0 

0. 003720 

0.995600 

1 .60 

0 .002  1  37 

0.002660 

0.998060 

L  .  50 

0  .  C  0  1 1  5  5 

0. CO  12 8  0 

0.999320 

1.60 

0 . 0  C  C  5  7 1 

0.  0005  20 

0.999860 

1  .  70 

0 .000266 

0.000160 

1 . 000000 

1 .80 

0.000083 

0. C000G0 

l .000000 

1.90 

0  .  CCOO  18 

0. CO 00 00 

l . occooo 

2.CC 

0 . 0  C  C  0  0 1 

0. 00 00 00 

1 .000000 

p 

d2 

Results  for 

Expected 

xi+l  E  513x 

Observed 

±  (mod  235-31) 

Accumulated 

d^ 

Frequency 

Frequency 

Frequency 

C.  1C 

0 .236832 

0.235160 

0.235160 

0 . 20 

0  .  1  76973 

0.  1  7  5080 

0.610260 

0.30 

C .  1  396  9  5 

0.  1  39880 

0.950120 

0.60 

C.  1  12718 

0.  1  10960 

0 . 66  1080 

0.50 

0. 090969 

0. 08 76 00 

0.  76  8680 

0.60 

0.072616 

0.0  76600 

0.823080 

0 . 70 

C .056769 

0.056920 

0.680000 

0.80 

0.062813 

0. 066960 

0.926960 

0.9C 

0.030631 

0.031120 

0.956080 

i.GC 

0.019333 

0. 019660 

0.975520 

1  .  1C 

o’ .  0  1  0  /  7  7 

0. 010660 

0.985960 

1. 20 

0.CC6365 

0.005920 

0.991880 

1 . 30 

0.003760 

0. 003760 

0.995660 

1 . 60 

0  .  C  0  2  1  3  7 

0.002280 

0.997920 

1 . 50 

0.001155 

0 .  C  0  1  2  0  0 

0.999120 

1.60 

0  .  C  C  0  5  7  1 

0. 0C0660 

0.999760 

1 .70 

0.CG0266 

0. CO 02 00 

0.999950 

1 . 80 

0.000083 

0.000060 

1 .000000 

1 .90 

0. COCO  18 

0. 0000 00 

1  .  OOCOOO 

2.CC 

C.CCCCOl 

0. 000000 

1 . OOCOOO 

d2 

Results  for 

x  =  513x 

Xi+1  "  ^ 

^  (mod  233) 

(mod  233) 
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Juncosa's  Normal  Conversion  for  16  Uniform  Deviates 
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Juncosa's  Normal  Conversion  for  6^  Uniform  Deviates 
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Marsaglia’s  Normal  Conversion  -  Part  1 
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Marsaglia’s  Normal  Conversion  -  Part  2 
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GAME  NUMBER 

* 

NUMBER  CF  TRIALS 

* 

GROSS  WINNINGS 

X 

or 

x  ^ 

j{:  $  $  $  £  3};  3{j  $  #  s}t 

X  X  -X. 
O'  O'  O' 

>S  X  X  X  X  X  X  X  J.  X  \*»  X  X  s',  J,  X  X  J. 

or  o'  o-  o-  V  o*  O'  v  o*  V  o'  O'  o*  V  v  o'  o'  o*  o 

-  X  X 
'  or  o' 

❖ 

l 

* 

66 

X 

r 

35 

2 

$ 

5 

£ 

0 

❖ 

❖ 

3 

* 

5 

* 

0 

* 

r 

4 

❖ 

9 

X 

or 

2 

❖ 

•># 

or 

5 

* 

14 

X 

r 

5 

❖ 

❖ 

6 

X 

or 

236 

O' 

126 

* 

* 

7 

* 

76 

$ 

3  7 

❖ 

$ 

a 

X 

or 

7 

❖ 

1 

❖ 

x 

or 

9 

* 

64 

X 

or 

31 

❖ 

* 

10 

* 

23 

X 

or 

10 

❖ 

or 

1  i 

* 

5 

❖ 

0 

j, 

V 

12 

❖ 

9 

O' 

2 

❖ 

V 

13 

V 

10 

* 

3 

X 

O' 

r 

14 

$ 

3  85 

X 

or 

208 

❖ 

X 

or 

15 

* 

18 

❖ 

7 

❖ 

v. 

V 

16 

* 

7 

❖ 

l 

X 

O' 

V 

l  I 

* 

60 

❖ 

32 

❖ 

j, 

or 

18 

* 

9 

$ 

2 

* 

« 

19 

£ 

28 

* 

13 

X 

O' 

V 

20 

❖ 

42 

❖ 

21 

❖ 

or 

21 

* 

15 

r 

5 

❖ 

❖ 

22 

* 

165 

❖ 

88 

❖ 

* 

23 

* 

16 

❖ 

6 

❖ 

j, 

or 

2  4 

& 

37 

❖ 

1  7 

❖ 

X 

O'- 

25 

* 

5 

❖ 

0 

❖ 

X 

o' 

26 

* 

30 

❖ 

13 

❖ 

X 

r 

27 

-X 

or 

7 

* 

1 

❖ 

X 

O' 

28 

X 

or 

27 

❖ 

12 

❖ 

£ 

29 

* 

9 

X- 

O' 

2 

X 

V 

❖ 

30 

* 

19 

❖ 

8 

X 

or 

31 

24 

❖ 

10 

❖ 

❖ 

32 

5 

❖ 

0 

X 

O' 

X 

o' 

33 

* 

35 

X 

V 

16 

X 

O' 

❖ 

34 

* 

254 

* 

133 

* 

35 

❖ 

11 

❖ 

3 

❖ 

* 

36 

X 

17 

❖ 

6 

X 

O' 

❖ 

37 

* 

3  44 

❖ 

184 

❖ 

❖ 

38 

* 

18 

X 

O' 

7 

❖ 

❖ 

39 

$ 

24  3 

X 

O' 

127 

X 

O' 

O' 

40 

* 

56 

❖ 

27 

* 

O' 

41 

* 

21 

* 

8 

❖ 

* 

42 

* 

10 

o* 

3 

J, 

T 

❖ 

43 

£ 

*  74 

❖ 

37 

❖ 

❖ 

4  4 

❖ 

35 

❖ 

L  5 

❖ 

X 

or 

45 

A 

55 

X 

O' 

2  7 

4= 

❖ 

46 

❖ 

l  71 

❖ 

8  7 

❖ 

X 

or 

4  7 

22 

❖ 

9 

❖ 

* 

48 

* 

13 

X 

O' 

4 

X 

or 

❖ 

4  9 

* 

23 

❖ 

9 

* 

X 

or 

50 

❖ 

82 

❖ 

41 

O' 

3jc 

*  £  *  4  4  %  $  $  ❖  *  £  *  £ 

£  y,:  3}  3ji  j|:  j|:  3^  3,;  •$  Jj:  jJ;  $  jJ:  ){;  sfc  jJ;  3|:  3|c 

V  O' 

>r  #  *  #  *  ❖  ❖  V  £  A:  ❖  v  #  #  £ 

<5*# 

AVERAGE  NUMBER 

CF 

TRIALS  UNTIL  RUIN 

IS 

58.42 

Chuck-a-Luck  Results  for  $5 
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vV 

"I* 

GAME  NUMBER 

•X 

't 

NUMBER  C  F  TRIALS  * 

GROSS  WINNINGS 

if 

ar 

❖  if  $  v  ^  i  ^  ^  $  -Jy  *  i  #  £  v  -r  ^ 

^  ^  Y  Y  Y  Y  Y  ^'v  Y  Y  ^  V  V  &  ^  ^  -fa  ^ 

if  if  5|o;;  if  3|;  if  if  >f  if  if  if 

if  if  if 

❖ 

l 

j. 

26  * 

8 

❖ 

2 

« 

39  * 

15 

JU 

Y 

❖ 

3 

* 

2  98  * 

157 

❖ 

❖ 

A 

$ 

35  * 

13 

Y 

* 

5 

a. 

T 

24  * 

7 

a. 

Y 

❖ 

6 

* 

73  * 

33 

* 

J, 

•V 

7 

* 

78  * 

36 

if 

i* 

8 

❖ 

74  * 

35 

a, 

Y 

❖ 

9 

* 

117  * 

61 

Y 

¥ 

10 

❖ 

225  * 

117 

if 

❖ 

1 1 

* 

183  * 

93 

>** 

'V 

❖ 

12 

* 

25  * 

8 

❖ 

❖ 

13 

if 

94  * 

46 

if 

❖ 

14 

if 

60  * 

27 

if 

❖ 

15 

'T 

18  * 

4 

if 

❖ 

16 

if 

16  * 

3 

if 

* 

17 

£ 

v  30  * 

1 1 

if 

£ 

18 

* 

44  * 

19 

if 

❖ 

19 

a. 

r 

128  * 

63 

if 

❖ 

20 

if 

122  * 

60 

ff 

❖ 

21 

if 

4  35  * 

230 

if 

❖ 

22 

if 

424  * 

225 

❖ 

V 

23 

❖ 

70  * 

32 

X 

■f 

•JU 

'¥• 

24 

❖ 

89  * 

43 

V 

v»„ 

V 

25 

* 

5  7  * 

26 

❖ 

J- 

V 

26 

if 

127  * 

63 

27 

* 

3  7  * 

14 

*f 

❖ 

28 

* 

134  * 

65 

if 

2  9 

* 

104  * 

52 

3- 

v 

❖ 

30 

* 

20  * 

5 

if 

3  l 

a* 

126  * 

63 

if 

« 

32 

❖ 

92  * 

43 

if 

❖ 

33 

* 

26  * 

8 

X 

V 

34 

* 

158  * 

80 

❖ 

❖ 

35 

* 

12  * 

1 

❖ 

❖ 

36 

2  31  * 

124 

❖ 

❖ 

37 

* 

2  79  * 

146 

*v 

¥ 

o. 

38 

* 

4  6  * 

19 

*> 

r 

❖ 

39 

186  * 

95 

❖ 

❖ 

40 

s 

56  * 

25 

❖ 

❖ 

41 

209  * 

109 

❖ 

❖ 

42 

* 

184  * 

95 

❖ 

❖ 

43 

=4= 

127  * 

63 

❖ 

❖ 

44 

* 

157  * 

80 

❖ 

❖ 

45 

❖ 

61  * 

28 

❖ 

* 

46 

if 

76  * 

3  7 

if 

47 

if 

15  * 

3 

a* 

Y 

❖ 

48 

if 

49  * 

20 

❖ 

❖ 

49 

if 

465  * 

2  47 

❖ 

if 

50 

if 

232  * 

l  19 

❖ 

y,:  ^5  if  ^  #  £  if  £  ^  ^  if  $ 

V  v  ^  ^  V  V  Y  Y  Y  Y  Y  ^  'r  V  Y  3!'  V  Y  V 

if  i|c  if  if  if  if  y:  if  if  if  : 

%«-  JU 
T  Y  Y 

AVERAGE  NUMBER 

CF 

TRIALS  UNTIL  RUIN  IS 

119*86 

Chuck-a-Luck  Results  for  $10 
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*  GAME  NUMBER 

vi, 

Y 

NUMBER  CF  TRIALS 

si. 

GROSS  WINNINGS 

>v 

Y 

^  ^  ^  ^  Y  ^  ^  -V  V  T  ^ 

-JC  J. 

Y  Y 

J.  Jj  J.  si  .*.  s«.  S.V  -A,  -A.  ,C  S^  O'  ^  f  si  V'  ^  ^  ^ 
Y  Y  Y  Y  Y  Y'  Y  Y  Y  Y  Y  Y  Y  Y  Y  Y  Y  Y'  V 

-s*.  A  x». 

Y  Y  *Y 

*4.  v1.  v*-  4/  4*  y<  4r  4»  4.  J.  4f  «l,  4,  «4^»  4/ 

Y  'i-  -V  'i'  Y  V  V  V  *i'  Y  Y  Y  vr  Y  Y 

*  1 

$ 

138 

63 

❖  2 

❖ 

430 

❖ 

220 

4. 

Y 

*  3 

❖ 

126 

>4, 

Y 

58 

4* 

Y 

£  4 

* 

96 

J. 

Y 

4  1 

4, 

Y 

*  5 

* 

83 

sV 

Y 

34 

4U 

Y 

#  6 

* 

267 

>4U 

Y 

137 

❖ 

❖  7 

$ 

2057 

■sV 

Y 

1107 

Y 

❖  8 

4. 

Y 

123 

s^ 

Y 

55 

v«, 

Y 

❖  9 

* 

114 

-I** 

51 

4. 

Y 

*  10 

* 

94 

Jt, 

Y 

4  0 

4. 

Y 

*  1 1 

* 

4  16 

sH> 

Y 

2  1  1 

nV 

Y 

*  12 

* 

132 

V 

60 

❖ 

*  13 

* 

356 

❖ 

187 

-A- 

'1* 

*  14 

$ 

217 

105 

❖ 

❖  15 

* 

5  43 

■v* 

289 

4, 

Y 

*  16 

« 

145 

* 

68 

4< 

Y 

*  17 

* 

466 

slU 

2  43 

❖ 

*  18 

❖ 

7  5 

Y 

3  l 

❖ 

*  19 

Y 

208 

vl, 

Y 

99 

4/ 

Y 

*  2  0 

* 

4  2 

❖ 

12 

4. 

Y 

*  21 

N 

64 

Y 

24 

Y 

*  22 

vC 

Y 

181 

❖ 

86 

Y 

*  23 

* 
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* 

137 

*  2  4 

$ 

195 

sl» 

'I* 

95 

❖ 

❖  25 
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❖ 

76 

4. 

Y 
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$ 
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V 
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Y 
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❖ 
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❖ 
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4. 

Y 
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❖ 
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4/ 

“V 
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❖ 
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* 
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❖ 

66 

V 
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-4 
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❖ 
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4r 

Y 
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V 
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❖ 

49 

V 
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* 
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❖ 

58 

4/ 

Y 
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* 
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❖ 

43 

A. 
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*  34 

❖ 
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❖ 

4  7 

❖ 
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❖ 
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Y 
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Y 

£  3  6 
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Y 

60 

*  37 

24  8 

❖ 

123 

❖ 

*  3  8 

* 
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V 

73 

❖ 

*  39 

❖ 

2  40 

4. 

Y 

1  19 

-.O 

Y 

*  4  0 

❖ 
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❖ 

79 

* 

❖  4  l 

❖ 

68 

^1' 
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❖ 

*  42 

* 
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4. 

Y 
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* 
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❖ 
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Y 
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❖ 
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❖ 

34 

❖ 
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CF 

TRIALS  UNTIL  RUIN 

I  s 

240.06 

Chuck- a- Luck  Results  for  $20 


• 

* 

»■ 

* 

108 


Throws 

Crosses 

1CCC0 

3334 

2CC00 

6  6  5  1 

30000 

9980 

4  C  C  C  C 

1  32  72 

5C0C0 

16595 

6  COCO 

19897 

7C00C 

2  319  3 

8  COCO 

26552 

9  C  CO  0 

29836 

1CCCC0 

33  176 

Result  s 

„  r  1  2 

for  J  x- 

0 

Throws 

Crosses 

1C00C 

2463 

2C0C0 

4956 

3CCC0 

74  5  2 

4CCC0 

9  9  3  3 

5C0GC 

12420 

6C0CC 

14  8  18 

7C0C0 

1  7  2  7  9 

8C0C0 

19  8  4  5 

9CC00 

22  3  10 

1CCCCC 

2  4  8  06 

Results  for 

0 


Area 

C  .  3334CC00 
0*332  55000 
0.33266667 
0*  331  8  0 COO 
0 .33190000 
0.33161667 
0.331  3285  7 
0. 331 9CC00 
0.33151111 
0.33176000 


dx  Using  Buffon’ 


Area 

0. 2463C000 
0.24780000 
0 . 2484C0Q0 
0.24832500 
0.24840 COO 
0.24696667 
0.24684286 
0.24806250 
0.24788889 
0. 248C6000 


dx  Using  Buffon’ 


Technique 


Technique 
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